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19.  (Concluded) 

DBGADIS  end  Gaussian  pluna  nodal  pradletlons.  (DEGADIS  is  an  atisospharic 
dispersion  model  designed  to  account  for  the  influences  of  denser-than-air 
gases.)  In  addition,  analysis  of  the  experimental  data  indicated  heat 
transfer  to  the  aerosol  cloud  was  ins ign^ leant  although  the  cloud  tem¬ 
perature  was  as  low  as  -60  centigrade.  mer 

Several  simplifying  assumptions  were  tried  when  determining  the  aerosol 
density  as  a  function  of  cloud  concentration  including:  (1)  a  pseudobinary 
mixture  of  ammonia  and  air/water  which  ignored  water  phase  changes ; 

(2)  ideal  ammonia/water  liquid  solution  behavior;  and  (3)  assumption  of  a 
linear  relationship  between  ammonia  concentration  and  mixture  density. 

Based  on  calculations  with  DEGADIS,  the  assumption  of  a  pseudobinary 
mixture  of  ammonia  and  air/water  proved  to  be  inadequate.  The  DEGADIS - 
predicted  distance  to  given  concentration  levels  were  calculated  asstimlng 
ideal  solutions  in  the  liquid  phase  for  the  ammonia/water  system;  these 
computations  reproduced  the  DEGADIS -predicted  distances  using  the  unsim¬ 
plified  thermodynamics  treatment.  Finally,  the  DEGADIS -predicted  distance 
to  given  concentration  levels  were  calculated  assuming  a  linear 
relationship  between  mixture  density  and  contaminant  concentration;  these 
computations  reproduced  the  DEGADIS -predicted  distances  using  the 
unslmpllfled  thermodynamics  treatment  to  within  about  20  percent.  In  the 
absence  of  heat  transfer,  this  last  assumption  may  be  useful  as  a  screening 
technique  for  hazard  assessment. 
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SECTION  I 
imODUCTIOH 

A.  OBJECTIVE 

The  objective  of  this  project  was  to  calculate  the  release  Richardson 
number  for  four  field-scale  releases  of  liquid  anhydrous  ammonia  to  determine 
the  Importance  of  cloud  density  In  the  atmospheric  dispersion  of  toxic  clouds 
resulting  from  these  test  spills.  After  calculating  the  release  Richardson 
number,  the  OEGAOIS  model  was  used  to  slanilate  the  atmospheric  dispersion  of 
the  toxic  clouds.  The  model  results  were  compared  to  field  measurements  to 
determine  the  suitability  of  DEGAOIS  to  simulate  these  types  of  releases. 

B.  BACKGROUHD 

The  U.  S.  Air  Force  (USAF)  has  been  taking  steps  to  Improve  the  safety 
procedures  available  during  operations  Involving  hazardous  chemicals  and 
fuels.  In  support  of  this  effort,  the  Air  Force  Engineering  and  Services 
Laboratory  (AFESC/RDV)  has  sponsored  research  on  the  dispersion  of 
denser-than-alr  gas  clouds  resulting  from  accidental  releases  of  toxic 
chemicals.  This  study  examines  four  field-scale  releases  of  liquid  anhydrous 
ammonia  from  the  "Desert  Tortoise"  (DT)  series  (Reference  1). 

Denser-than-alr  gases  are  of  particular  Interest  because  current  (passive) 
dispersion  models  used  to  predict  concentrations  (and  safety  limits)  may  be 
Inappropriate.  Models  currently  used  for  safety  assessment  assume  that  the 
dispersing  contaminant  does  not  affect  the  atmospheric  flow  field.  This  may 
not  be  true  for  denser-than-alr  gases;  dispersion  rates  may  be  considerably 
reduced  for  denser-than-alr  gases  due  to  the  stable  stratification  which  may 
be  present.  The  release  Richardson  ntaiber  has  been  used  successfully  as  a 
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criterion  for  determining  whether  a  particular  release  of  material  can  be 
adequately  modeled  using  passive  dispersion  techniques  or  whether  denser-than- 
alr  gas  effects  snist  be  considered.  (The  release  Richardson  number  is  a 
dimensionless  number  relating  the  effects  of  the  type  of  material  released, 
the  ambient  wind  profile,  and  characteristics  of  the  release  Itself.)  This 
study  examines  the  applicability  of  the  release  Richardson  number  to  the  DT 
tests.  The  liquid  anhydrous  ammonia  was  released  as  a  pressurized  liquid  Jet 
in  the  DT  tests;  this  study  also  attempts  to  quantify  the  effect  of  the  Jet 
behavior  on  the  relative  Importance  of  the  denser-than-alr  gam  effects  present 
In  the  releases. 

C.  SCOPB/APPROACH 

Because  the  liquid  anhydrous  ammonia  was  released  as  a  pressurized  Jet, 
the  ammonia  formed  a  cold  ammonia  aerosol  as  it  moved  downwind.  Upon  mixing 
It  with  the  ambient  air,  some  of  the  ambient  humidity  condensed  to  form  a 
nonideal  liquid  phase  composed  of  ammonia  and  water.  This  study  examines  the 
influence  on  the  dispersion  processes  by  the  aerosol  formation,  the  nonideal 
liquid  solution  behavior  of  the  ammonia  and  water,  and  the  heat  transfer  to 
the  cold  aerosol;  potential  simplifications  are  discussed. 

Finally,  this  study  examines  the  consistency  of  model  predictions  using 
the  Gaussian  plume  model  and  DB6ADIS  (a  computer  model  designed  to  account  for 
the  influences  of  denser-than-alr  gases)  with  the  results  of  the  DT  tests. 

The  importance  of  such  factors  as  the  aerosol  behavior,  heat  transfer  to  the 
aimnonla  cloud,  and  initial  Jet  behavior  on  the  model  predictions  are 
discussed.  The  sensitivity  of  DEGADIS  to  uncertainty  in  input  parameters  not 
routinely  available  at  operational  USAF  sites  is  also  assessed. 
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SECTION  II 
METHODOLOGY 

Four  field-scale  releases  of  liquid  anhydrous  ammonia  (the  "Desert 
Tortoise”  (DT)  series)  were  performed  in  1983  by  La%rrence  Livermore 
National  Laboratories  (LLNL)  for  the  O.S.  Coast  Guard,  the  Fertilizer 
Institute,  and  Environment  Canada  (Reference  1).  The  releases  were 
conducted  on  Frenchman  Flat,  a  dry  lake  bed  approximately  6  kilometers  long 
and  3  kilometers  vide,  at  the  Nevada  Test  Site.  Although  the  area  is 
normally  very  dry,  unusually  heavy  rains  had  occurred  before  the  test 
period.  Water  was  standing  on  the  test  site  during  the  first  three  tests, 
but  the  lake  bed  was  dry  during  Test  4.  A  summary  of  the  release 
conditions  for  the  tests  is  presented  in  Table  1.  A  surface  roughness  of 
0.003  meters  was  used  for  the  site  (Reference  1).  Pressurized  storage 
resulted  in  rapid  expansion  of  the  ammonia  as  it  exited  through  an  orifice 
at  the  end  of  the  spill  pipe .  This  section  summarizes  important  aspects  of 
this  type  of  release  specifically  dealing  with  the  initial  properties  of 
the  ammonia  aerosol  which  was  formed.  Also,  the  dispersion  processes 
downwind  of  the  source  are  discussed.  The  succeeding  sections  address  the 
analysis  of  the  DT  test  data  to  yield  information  suitable  for  model 


TABLE  1.  SUMMARY  OF  RELEASE 

CONDITIONS  FOR 

THE  DESERT 

TORTOISE 

SERIES 

DTI 

DT2 

DT3 

DT4 

Spill  rate  (kg/s) 

81. 

117. 

133. 

108. 

Windspeed  at  5.83  m  (m/s) 

8.33 

6.01 

8.16 

4.99 

Friction  velocity  (m/s) 

0.442 

0.339 

0.448 

0.286 

Ambient  temperature  (”C) 

28.7 

30.9 

33.8 

33.2 

Ambient  pressure  (bar) 

0.909 

0.910 

0.907 

0.903 

Relative  humidity  (%) 

13.2 

17.5 

14.8 

21.3 

Monin-Obukhov  length  (m) 

92.7 

94.7 

570.7 

45.2 

Fasquill  stability  class 

D 

D 

D 

E 
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comparison,  a  method  of  determining  the  relative  importance  of  the  various 
flow  and  dispersion  processes  present,  and  the  comparison  of  the  results 
from  the  DT  test  data  with  a  Gaussian  plume  model  and  DEGADIS. 

A.  AMMONIA  AEROSOL  FORMATION  AND  SUBSEQUENT  MIXING  WITH  AIR 

The  release  of  anhydrous  ammonia  from  storage  at  elevated  pressure  and 
ambient  temperature  results  in  the  formation  of  a  denser- than -air  ammonia 
aerosol  if  the  release  is  from  the  liquid  phase;  such  releases  are 
typically  violent  jets.  The  DT  release  conditions  represent  this  type  of 
release.  (A  general  discussion  of  the  behavior  of  the  early  phases  of 
ammonia  releases  is  presented  by  Wheatley  (References  2  and  3);  for  this 
report,  only  the  conditions  of  the  DT  releases  are  discussed.)  The 
behavior  of  the  flashing  liquid  jet  can  be  described  (approximately)  as 
either  isenthalpic  or  isentropic.  (Isenthalpic  behavior  implies  that  all 
of  the  energy  of  the  liquid  jet  goes  into  the  energy  of  the  resulting 
aerosol,  ignoring  any  kinetic  energy;  isentropic  behavior  implies  that  the 
flow  is  adiabatic,  frictionless,  and  ignores  other  irreversibilities  such 
as  shocks.)  For  storage  temperatures  below  about  SO’C,  the  mass  fraction 
of  liquid  ammonia  in  the  resulting  aerosol  is  approximately  the  same 
regardless  of  whether  the  flow  is  assumed  to  be  isenthalpic  or  isentropic. 
Assuming  isenthalpic  flow,  the  aerosol  was  estimated  to  have  (initially)  a 
temperature  of  -36°C  with  about  81  percent  (by  mass)  of  ammonia  in  the 
liquid  phase.  For  hazard  assessment  purposes,  the  distance  required  for 
the  flashing  liquid  jet  to  regain  thermal  equilibrium  and  form  the  aerosol 
is  insignificant  (on  the  order  of  1  meter) . 

Since  dispersion  depends  on  the  density  of  the  air/ammonia  mixture,  it 
is  necessary  to  estimate  the  mixture  density  as  a  function  of  ammonia  mole 
fraction.  Some  of  the  ambient  humidity  will  condense  when  air  is  entrained 
into  the  aerosol.  Since  ammonia  and  water  form  nonideal  solutions  in  the 
liquid  phase,  consideration  of  the  nonideal  behavior  of  the  liquid  phase 
may  be  important  under  some  circumstances.  For  the  ammonia/water  system, 
Wheatley  (References  4  and  5)  has  investigated  these  interactions  and 
incorporated  them  into  a  computer  code,  TRAUMA,  which  was  used  to  estimate 
the  ammonia/air/water  mixture  density  as  a  function  of  ammonia  concentra¬ 
tion  for  the  DT  releases.  Figure  1  shows  the  mixture  density  and  mixture 
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Ammonia  Mole  Fraction  Ammonia  Mole  Fraction 

Figure  1.  Ammonle/Air/Vater  Mixture  Temperature  and  Density  as 
Ftinctlons  of  Ammonia  Mole  Fraction  for  DT4  Conditions 
Calculated  using  TRAUMA. 


temperature  as  functions  of  ammonia  mole  fraction  for  DT4  which  Is  typical 
of  all  of  the  releases  calculated  using  TRAUMA  version  1.0.  (The  equilib¬ 
rium  mixture  temperature  decreases  with  the  addition  of  air.) 

Three  simplifying  assumptions  can  be  made  to  describe  the  ammonla/alr/ 
water  mixture  density:  (1)  a  pseudobinary  mixture  of  ammonia  and  air/ 
water;  (2)  Ideal  ammonia/water  liquid  solution  behavior;  and  (3)  assumption 
of  a  linear  relationship  between  ammonia  concentration  (In  kg/m*')  and  mix¬ 
ture  density.  To  determine  whether  these  assumptions  are  appropriate,  the 
effect  on  DEGADIS  predictions  of  the  DT4  release  conditions  was  examined. 
TRAUMA  was  used  to  determine  the  mixture  density  of  a  pseudobinary  mixture 
of  alr/water  and  ammonia;  water  phase  changes  were  Ignored,  greatly  sim¬ 
plifying  the  calculation  procedure  required  (since  the  partial  pressure  of 
ammonia  Is  then  eq\ial  to  Its  vapor  pressure) .  Ammonia  phase  changes  were 
taken  Into  account.  (For  details  of  the  capabilities  of  TRAUMA,  see 
Wheatley  (Reference  6).)  Comparing  simulations  made  with  this  simplified 
density  specification  and  the  density  specification  of  Figure  1,  substan¬ 
tial  differences  were  noted  In  the  DEGADIS  dispersion  calculation  for 
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ammonia  concentrations  below  on  the  order  of  1  percent  for  the  release 
conditions  of  DT4.  TRAUMA  was  also  used  to  determine  the  ammonla/alr /water 
mixture  density  assuming  Ideal  liquid  solution  behavior;  this  assumption 
also  simplifies  the  calculation  procedure  required  (since  the  partial 
presstire  of  ammonia  is  then  equal  to  its  vapor  pressure  times  its  liquid 
phase  mole  fraction) .  Heat  of  mixing  effects  between  the  ammonia  and  water 
in  the  liquid  phase  were  ignored.  The  assumption  of  ideal  liquid-phase 
solution  behavior  adequately  described  the  mixture  density  so  that  no 
difference  was  noted  in  the  DEGADIS  dispersion  calculations  based  on  the 
release  conditions  of  DT4  (for  relative  humidities  ranging  from  0  to  100 
percent).  Finally,  a  linear  relationship  between  ammonia  concentration  (in 
kg/m^)  and  mixture  density  was  used  directly  in  DEGADIS;  the  endpoints  were 
represented  by  the  released  ammonia  aerosol  and  the  ambient  air.  Upon 
comparison,  the  DEGAOIS-predicted  distances  using  this  simplified  scheme 
were  within  about  +20  percent  of  the  DEGADIS-predicted  distance  to  a  given 
concentration  level  using  the  density  specification  of  Figure  1;  this 
agreement  was  checked  to  about  100  ppm  for  the  conditions  of  DT4. 

B.  JET  AMD  NONJET  FLOW  PHASES 

The  initial  growth  of  a  turbulent  jet  is  characterised  by  large-scale 
turbulent  motions  which  entrap  air  in  the  Jet  flow  (Reference  7).  This 
initial  near-jet  region  was  estimated  to  persist  on  the  order  of  10  meters 
for  the  DT  series.  As  a  result  of  air  entrainment  and  expansion  of  the 
jet,  the  momentum  of  the  jet  decreases,  and  the  jet  is  well  described  by 
self-similar  velocity  and  concentration  profiles  in  the  absence  of  inter¬ 
action  with  the  ground  (Reference  8).  The  jet  velocity  decreases  (due  to 
air  entrainment)  and  approaches  the  ambient  wind  velocity.  For  the  DT 
series,  the  ammonia  was  released  0.79  meters  above  ground  level,  and  the 
jet  struck  the  ground  within  one  meter  or  so.  Therefore,  the  DT  releases 
would  be  expected  to  behave  like  a  ground-level  (jet)  release  due  to  the 
interaction  with  the  solid  boundary.  In  contrast  to  nonaerosol  jets,  the 
aerosol  jets  in  the  DT  releases  may  be  subject  to  rainout  (deposition  of 
some  of  the  liquid  phase  on  the  ground).  Wheatley  (Reference  5)  has 
examined  the  possibilities  of  rainout  for  ammonia  aerosols  and  found  that 
rainout  is  not  expected  for  pressurized  releases  in  the  absence  of  solid 
boundary  interactions.  For  the  DT  series.  Goldwire  et  al.  (Reference  1) 
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reported  thet  some  of  the  ammonie  was  deposited  on  the  groiind  and  formed  a 
liquid  pool,  probably  because  the  release  vas  near  groxind  level. 

After  the  contaminant /air  mixture  Jet  velocity  has  diminished,  non jet 
flow  phases  become  Important  Including:  (1)  negative  buoyancy -dominated 
flow;  (2)  stably  stratified  shear  flow;  and  (3)  passive  dispersion  due  to 
atmospheric  flow.  The  theory  underlying  the  dispersion  prediction  of  trace 
contaminants  (passive  dispersion)  generally  assumes  that  the  dispersion  is 
the  result  of  atmospheric  turbulence.  Although  characterization  of  the 
atmospheric  flow  suffers  from  the  limits  of  understanding  of  turbulent 
fluid  motion,  there  Is  a  fairly  well-developed  theoretical  basis  for  pre¬ 
diction  of  passive  atmospheric  dispersion,  along  with  an  extensive  experi¬ 
mental  data  base  derived  from  atmospheric  flow  measurements  (Reference  9) . 
In  contrast  to  passive  dispersion,  the  release  of  large  quantities  of  a 
denser-than-alr  gas  (DTAG)  Into  the  atmosphere  can  significantly  alter  the 
atmospheric  flow  In  the  vicinity  of  the  release.  For  some  releases,  this 
negative  buoyancy -dominated  (gravity-driven)  flow  and  the  resulting 
Interaction  with  the  atmospheric  flow  can  have  an  Important  effect  on  the 
distance  needed  to  reduce  the  concentration  to  a  given  level.  Between 
these  two  extremes,  the  stably  stratified  shear  flow  phase  differs  from 
neutrally  buoyant  flow  and  Is  characterized  by  the  following:  (1)  a 
crosswind  gravity -driven  flow  due  to  the  negative  buoyancy  of  the  flow  Is 
present;  (2)  because  of  this  gravity-driven  spreading,  these  pliimes  tend  to 
be  wider  and  shallower  than  a  neutrally  buoyant  plume  under  the  same 
conditions;  and  (3)  because  of  stable  vertical  density  stratification,  the 
vertical  mixing  rate  Is  reduced.  (A  detailed  discussion  of  each  of  these 
nonjet  flow  phases  Is  found  In  Havens  and  Spicer  (Reference  10);  see 
Section  III  and  Appendix  A  for  a  method  to  estimate  the  relative  Importance 
of  the  jet  and  nonjet  flow  phases.) 
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SECTION  III 

ANALYSIS  OF  DESERT  TORTOISE  DATA  FOR  MODEL  COMPARISON 


For  the  DT  releases,  two  sensor  arrays  were  used  to  measure  the  down- 
vind  ammonia  concentrations.  The  mass  fl\ix  array  consisted  of  a  row  of 
seven  gas  sensor  masts  located  100  meters  downwind  of  the  release  area. 

The  primary  gas  concentration  sensors  in  this  array  were  Mine  Safety 
Appliance  (MSA)  nondispersive  IR  gas  sensors  which  were  heated  to  vaporize 
the  aerosol  so  that  the  total  ammonia  concentration  could  be  determined. 
(The  velocity  of  gas  through  the  100-meter  array  was  to  be  measured  with 
anemometers  to  determine  the  mass  flux  of  gas  through  the  array,  but 
corrosion  caused  by  the  ammonia  made  these  instruments  unreliable.)  At  800 
meters  downwind  of  the  release  area,  the  instrument  array  consisted  of  five 
sensor  masts  spaced  100  meters  apart;  each  mast  had  three  gas  sensors  and 
three  thermocouples.  The  sensors  were  located  at  1.0,  3.5,  and  8.5  meter 
heights  on  each  sensor  mast.  The  primary  gas  sensors  in  this  array  were 
International  Sensor  Technology  (1ST)  solid-state  gas  sensors.  In  addition 
to  these  sensor  arrays,  eight  portable  1ST  gas  sensors  were  located  at 
1 -meter  elevations  from  1.4  to  5.5  kilometers  downwind. 

To  facilitate  the  comparison  of  model  simulations  and  reported  concen¬ 
tration  and  temperature  measurements,  the  reported  concentrations  and 
temperatures  along  with  reported  wind  trajectories  and  mass  rates  at  800 
meters  were  analyzed  to  determine  the  best  representation  of  steady-state 
concentrations  and  temperatures.  (The  reported  wind  trajectories  did  not 
correspond  exactly  with  the  plume  centerline  due  to  the  momentum  of  the 
release.)  The  mass  rate  passing  the  800-meter  array  was  reported  by  LLNL 
on  the  basis  of  the  observed  concentration  and  temperature  profiles  at  the 
800-meter  array  and  velocity  profiles  Interpolated  from  other  meteorologi¬ 
cal  towers.  The  mass  rate  passing  the  800-meter  array  for  each  of  the 
tests  during  the  time  period  when  the  concentration  and  temperature  pro¬ 
files  were  examined  are  summarized  in  Table  2.  Most  of  the  released  mass 
is  accounted  for  at  the  800-meter  array  (67  percent  to  89  percent)  for  all 
tests  except  DT3  (42  percent) .  Ignoring  the  uncertainty  in  the  estimated 
ammonia  mass  rate  at  800  meters,  the  unaccounted  ammonia  mass  at  the  800- 
meter  array  could  have  been  deposited  on  the  ground  by  rainout  or  could 
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have  been  absorbed  by  the  water  on  the  ground  in  DTI  -  DT3.  Tables  3-6 
suimnarlze  the  reported  concentration  and  temperature  for  the  time 
reflecting  these  (approximate)  steady-state  conditions. 


TABLE  2.  REPORTED  AMMONIA  MASS  RATES  AT  800-METER 
SENSOR  ARRAY  AT  THE  SPECIFIED  TIMES 


Release 

Integrated  Contaminant 
Mass  Rate  at  800  m  (kg/s)* 

Time  from  Start 
of  Release  (s) 

DTI 

66.5 

-240. 

DT2 

104.6 

-270. 

DT3 

56.0 

-200. 

DT4 

72.1 

-340  to  390 

I  - 

★from  Reference  1. 


TABLE  3.  REPORTED  CONCENTRATION  AND  TEMPERATURE  FOR  DTI  AT  800  METERS 

AND  240  SECONDS  REFLECTING  (APPROXIMATE)  STEADY -STATE  CONDITIONS 


Sensor 

Elevation 

(m) 

Concentration 

(%)  /  Temperature  (°C) 

Sensor 

Station 

G24 

Sensor 

Station 

G23 

Sensor 

Station 

G22 

Sensor 

Station 

G21 

Sensor 

Station 

G20 

8.5 

0.0/29.1 

0.07/27.4 

0.02/28.6 

0.0/29.6 

0.0/29.6 

3.5 

0.0/28.9 

0.88/22.7 

0.88/22.7 

0.0/29.3 

0.0/29.3 

1.0 

0.0/28.8 

1.00/24.0 

1.00/24.0 

0.0/29.4 

0.0/29.4 

Note;  The  plume  centerline  was  approximately  midway  between  G22  and  G23. 
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TABLE  4.  REPORTED  CONCENTRATION  AND  TEMPERATURE  FOR  DT2  AT  800  METERS 

AND  270  SECONDS  REFLECTING  (APPROXIMATE)  STEADY- STATE  CONDITIONS 


Concentration 

(%)  /  Temperature  (“C) 

Sensor 

Elevation 

(m) 

Sensor 

Station 

G24 

Sensor 

Station 

G23 

Sensor 

Station 

G22 

Sensor 

Station 

G21 

Sensor 

Station 

G20 

8.5 

0 . 0/30 . 3 

0 . 14/28 . 0 

0.23/26.7 

0.0/30.4 

0.0/30. A 

3.5 

0.0/30.3 

1.08/25.7 

1.70/21.7 

0.91/23.3 

0.0/30.7 

1.0 

0.0/30.0 

1.67/25.0 

1.86/22.5 

1.71/22.2 

0.0/30.7 

Note:  The  plume  centerline  was  approximately  over  G22. 

TABLE  5. 

REPORTED  CONCENTRATION  AMD  TEMPERATURE  FOR  DT3  AT  800 
AND  200  SECONDS  REFLECTING  (APPROXIMATE)  STEADY-STATE 

METERS 

CONDITIONS 

Concentration 

(%)  /  Temperature  (®C) 

Sensor 

Elevation 

(m) 

Sensor 

Station 

G24 

Sensor 

Station 

G23 

Sensor 

Station 

G22 

Sensor 

Station 

G21 

Sensor 

Station 

G20 

8.5 

0.0/32.7 

* 

0.0/33.4 

0.0/33.3 

0.0/33.1 

3.5 

0.58/29. 

2  * 

0.78/29.2 

0.0/33.3 

0.0/33.3 

1.0 

0.87/28. 

1  * 

1.00/28.1 

0.0/32.9 

0.0/32.9 

Note:  The  plume  centerline  vas  approximately  over  G23. 
G23  failed  to  record  any  data  due  to  Instrument  failure. 
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TABLE  6.  REPORTED  CONCENTRATION  AND  TEMPERATURE  FOR  DT4  AT  800  METERS 

AND  350  SECONDS  REFLECTING  (APPROXIMATE)  STEADY-STATE  CONDITIONS 


Sensor 

Elevation 

(m) 

Concentration 

(%)  /  Temperature  (°C) 

Sensor 

Station 

G24 

Sensor 

Station 

G23 

Sensor 

Station 

G22 

Sensor 

Station 

G21 

Sensor 

Station 

G20 

8.5 

0.0/32.9 

0.0/32.0 

0.0/32.0 

0.0/32.2 

0.0/32.2 

3.5 

0.0/32.2 

0.77/25.9 

1.26/23.9 

0.80/24.8 

0.47/27.4 

1.0 

0.0/31.4 

1.29/26.3 

1.94/23.9 

1.90/23.3 

0.83/27.3 

Note:  The  plume  centerllni  was  approximately  midway  between  G21  and  G22. 


For  model  comparison  purposes,  the  measured  concentrations  were  fitted 
to  two  similarity  profiles: 


and 


(1) 


(2) 


where  c^  represents  the  maximum  centerline  contaminant  concentration,  z 
represents  the  height  above  ground,  y  represents  the  lateral  distance  from 
the  centerline,  and  are  Gaussian  (Pasquill-Gifford)  dispersion 
coefficients,  S^  and  S^  ere  distribution  coefficients  used  in  the  DEGADIS 
model,  and  a  Is  the  power-law  wind  profile  parameter.  The  concentration 
profiles  of  Equations  (1)  and  (2)  were  integrated,  along  with  a  (typical) 
velocity  profile  of  the  form 


la 


ll 


to  give  the  mess  rate  passing  a  plane  downwind: 

cu^dzdy  (A) 

Using  Equation  (4)  and  the  reported  mass  rates  passing  the  800-meter  array, 
ay,  Oj.,  Sy.  and  were  then  determined  by  a  weighted  nonlinear  least- 
squares  technique  from  the  reported  concentrations  and  temperatures  of 
Tables  3-6.  The  values  of  ay.  a^.  and  c^  for  Equation  (1)  and  Sy,  S^,  and 
Cj.  for  Equation  (2)  are  summarized  In  Tables  7  and  8,  respectively.  Notice 
that  the  maxlmiim  centerline  concentration  is  different  for  the  different 
profiles.  Furthermore,  notice  that  ay  =«  Sy/T?  and  a^  =»  S^  (to  within  about 
10  percent) . 


TABLE  7.  ESTIMATES  OF  Oy  AND 
a^  FOR  THE  DESERT 
TORTOISE  SERIES 
AT  800  METERS 


TABLE  8.  ESTIMATES  OF  S  AMO 
S^  FOR  THE  DESERT 
TORTOISE  SERIES 
AT  800  METERS 


a 

<m) 

^z 

<m) 

®c  3 

(kg/m^) 

DTI 

31.2 

3.42 

0.0267 

DT2 

65.7 

3.60 

0.0241 

DT3 

66.0 

2.50 

0.0155 

DT4 

86.7 

2.80 

0.0218 

Sy 

(m) 

Sz 

(m) 

®c  3 

(kg/m^) 

DTI 

43.8 

3.78 

0.0323 

DT2 

105.0 

4.80 

0.0205 

DT3 

106.2 

2.66 

0.0171 

DT4 

143.5 

3.11 

0.0222 
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SECTION  IV 

release  RICHARDSON  NUMBER 

The  purpose  of  this  section  Is  to  examine  one  way  of  determining 
whether  the  ground- level  jet  dominates  a  release  or  whether  the  other 
nonjet  phases  of  a  gro\ind-leveI  release  (including  the  negative  buoyancy- 
dominated  dispersion  phase,  the  stably  stratified  shear  flow  phase,  and  the 
passive  dispersion  phase)  are  the  only  important  dispersion  phases.  (The 
derivation  of  these  criteria  is  discussed  in  Appendix  A.) 

From  Havens  and  Spicer  (Reference  10)  and  Spicer  and  Havens  (Reference 
11),  the  previous  criteria  for  determining  which  phase  of  the  dispersion 
process  dominates  a  particular  ground-level,  low  initial  momentum  (nonjet) 
release  was  based  on  water  tunnel  experiments  reported  by  Britter  (Refer¬ 
ence  12).  In  that  sec  of  experiments,  a  salt/water  solution  was  released 
at  floor  level,  and  the  lateral  and  upwind  extent  of  the  plume- was 
recorded.  Analysis  of  Britter 's  releases  showed  the  following  criteria  of 
a  release  Richardson  ntimber  were  obtained  for  ground-level,  nonjet 
releases : 

If  Rig  5=  30  negative  buoyancy- dominated  phase 

If  1  ^  Rig  ^  30  stably  stratified  shear  flow  phase 

If  Rig  Z  1  passive  dispersion  phase 

where  Rig  -  g(PE  ‘  vertical  length  scale  H  is  approxi¬ 

mated  as  H  -  Q/uD  for  these  ground-level,  nonjet  releases.  (For  Britter 's 
releases,  the  ground- level  release  was  directed  upward  so  that,  in  the 
absence  of  entrainment  and  gravity  spreading  at  the  source ,  the  depth  of 
material  moving  downstream  (H)  would  just  be  Q/uO  for  a  uniform  approach 
velocity  u.)  Rig  represents  a  ratio  of  the  potential  energy  characteristic 
of  the  release  to  a  measure  of  the  ambient  turbulent  kinetic  energy.  For 
the  DT  series,  this  approximation  for  H  is  not  appropriate.  To  best 
approximate  the  potential  energy  characteristic  of  the  release  as  used  in 
Rig,  the  value  of  H  used  herein  for  the  DT  releases  is  2.0  meters  for  all 
the  releases  based  on  observation  of  the  photographic  records  of  the  tests 
as  well  as  the  height  of  the  release  pipe  outlet  (0.79  meters).  Using  this 
value  for  H,  the  values  of  Rig  for  the  DT  series  are  shown  in  Table  9. 

(For  estimation  of  Rig,  values  of  u*  were  taken  from  Goldwire  et  al. 
(Reference  1),  and  values  of  the  initial  aerosol  density  were  estimated 


TABLE  9.  VALUES  OP  Rt^  AND  (V/u*)^ 
FOR  THE  DT  SERIES 


Test 

Ric 

(V/u*) 2 

DTI 

280 

2300 

DT2 

490 

5000 

DT3 

270 

3700 

DT4 

630 

6100 

using  I«AUMA  as  discussed  In  Section  I.)  As  Indicated  by  Rl^.,  all  of  the 
DT  releases  would  be  In  the  negative  buoyancy- dominated  flow  phase  at  the 
release  In  the  absence  of  any  jet  effects. 

After  criteria  of  domination  of  a  particular  nonJet  dispersion  phase 
have  been  established,  It  can  be  determined  whether  a  particular  release  Is 
dominated  by  jet  effects.  Of  course,  jet  effects  would  dominate  a  release 
In  the  absence  of  wind  for  a  horizontal  release.  Goldwlre  et  al.  (Refer¬ 
ence  1)  suggest  that  jet  effects  were  Important  In  the  DT  series  to  over 
100  meters  downwind  of  the  release  because  agreement  between  the  position 
of  the  vapor  cloud  and  wind  trajectories  was  only  fair.  For  this  analysis, 
a  release  will  be  considered  dominated  by  jet  effects  associated  with  the 
release  If  the  rate  of  air  entrainment  due  to  the  jet  dominates  the  rate  of 
air  entrainment  due  to  the  dominant  nonjet  dispersion  phase  of  the  release. 
Based  on  the  analysis  In  Appendix  A.  ground-level,  horizontal  jet  effects 
dominate  the  dominant  non jet  flow  phase  when: 

V/u  5^  0.8  for  passive  dispersion  phase 

V/u  5  16/(19+Rig)  for  stably  stratified  shear  flow  phase 

(V/u*)^  5  10  Rig  for  negative  buoyancy- dominated  phase 

where  V  Is  the  velocity  of  the  material  leaving  the  pipe.  For  the  DT 
series,  the  appropriate  condition  to  check  Is  (V/u^)^  10  Rig  since  the 

dominant  nonjet  flow  phase  is  the  negative  buoyancy-dominated  flow  phase. 
Ratios  of  (V/u^)^  are  also  shown  In  Table  9.  Based  on  this  criteria,  the 
only  release  which  would  be  clearly  dominated  by  jet  effects  Is  DT3,  and 


the  reneinlng  tests 
properly  takes  Into 
above . 
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SECTION  V 

COMPARISON  OF  THE  PASQUILL-HANNA  GAUSSIAN  PLUME  MODEL 
AND  DEGADIS  WITH  THE  DESERT  TORTOISE  DATA 

The  Pasquill -Hanna  Gaussian  plvune  model  has  proven  to  be  applicable  to 
atmospheric  dispersion  problems  when  the  dispersion  of  the  contaminant  is 
only  a  function  of  the  atmospheric  turbulence  and  the  plume  does  not 
perturb  the  ambient  flow  field  (passive  dispersion).  The  steady-state 
Gaussian  plume  model  for  ground-level  releases  is  given  by 


where  y^  is  the  contaminant  mole  fraction,  Q  is  the  volumetric  source 
evolution  rate,  and  and  are  the  standard  deviation  of  the  lateral  and 
vertical  concentration  distributions,  respectively.  Values  of  were 
taken  from  Hanna  et  al.  (Reference  13),  and  values  for  172  which  take  into 
account  surface  roughness  effects  were  taken  from  Pasquill  and  Smith 
(Reference  9). 

The  DEGADIS  (DEnse  GAs  Dispersion)  model  was  developed  for  the  U.S. 
Coast  Guard  and  the  Gas  Research  Institute  and  was  designed  to  model  the 
atmospheric  dispersion  of  DTAG's  (References  10,  11).  DEGADIS  is  an 
adaptation  of  the  Shell  HEGADAS  model  described  by  Colenbrander  (Reference 
14)  and  Colenbrander  and  Puttock  (Reference  15);  DEGADIS  also  incorporates 
some  techniques  used  by  van  Ulden  (Reference  16) .  The  model  was  developed 
to  predict  the  dispersion  of  gas  from  a  ground- level  area  source  (such  as  a 
boiling  liquid  pool),  and  describes  three  phases  of  dispersion  which 
typically  occur  following  nonjet,  gro\md-level  DTAG  releases  discussed  in 
Section  II.  The  near-field,  buoyancy-dominated  dispersion  phase  is  modeled 
using  a  lumped  parameter  model  of  a  denser- than-air  "secondary  source" 
cloud  which  incorporates  air  entrainment  at  the  gravity- spreading  front 
using  a  frontal  entrainment  velocity;  this  description  is  based  on  labora¬ 
tory  data  reported  in  Havens  and  Spicer  (Reference  10) .  The  downwind 
dispersion  phase  of  the  calculation  assumes  a  power  law  concentration 
distribution  in  the  vertical  direction  and  a  modified  Gaussian  profile  in 
the  horizontal  direction  with  a  power  law  wind  profile  (Figure  2) .  The 
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vertical  mixing  rate  is  based  on  laboratory- scale  data  for  vertical  mixing 
in  stably  density-stratified  fluids  reported  by  Kantha  et  al.  (Reference 
17),  Lofquist  (Reference  18),  and  McQuaid  (Reference  19).  The  vertical 
dispersion  parameters  and  the  horizontal  dispersion  parameter  deter¬ 
mine  the  vertical  and  horizontal  profiles,  respectively.  It  should  be 
noted  that  the  rates  of  change  of  and  approach  the  rates  of  change  of 
cTy  and  a^,  respectively,  as  the  density  of  the  plume  approaches  the  ambient 
density.  A  complete  description  of  DEGADIS  is  included  in  Appendix  B;  a 
discussion  of  the  sensitivity  of  DEGADIS  to  uncertainty  in  model  input 
parameters  not  routinely  available  at  operational  USAF  sites  is  included  in 
Appendix  C. 

At  present,  DEGADIS  does  not  have  provisions  for  describing  jet 
releases.  However,  it  is  characteristic  of  DEGADIS  predictions  that  the 
distance  to  a  given  concentration  level  is  not  a  strong  function  of  the 
source  area  when  the  emission  flxix  is  greater  than  or  eqtial  to  the  maxinmm 
atmospheric  takeup  f Itix  (within  an  order  of  magnitude) .  Using  the 
arguments  presented  in  Section  IV,  the  diameter  of  the  source  used  in  the 
simulations  was  given  by  D  ■>  Q/uH  along  with  the  assumed  value  of  H  -  2 
meters  for  all  releases.  (Essentially  identical  values  for  c^,,  Sy,  S^,  and 
b  were  predicted  using  these  values  for  D  and  lOD.) 

In  addition,  the  thermodynamics  treatment  in  DEGADIS  does  not  presently 
provide  for  phase -change  heat  effects  (other  than  ambient  humidity  conden¬ 
sation)  in  the  gas/air  mixture.  Therefore,  aerosol  phase  change  effects  on 
the  cloud  density  were  accoimted  for  by  estimating  the  mixture  density  as  a 
function  of  contaminant  concentration  using  TRAUHA  as  discussed  in  Section 
II .  This  assumes  that  a  given  amount  of  pure  released  aerosol  and  ambient 
htimid  air  are  mixed  until  they  reach  equilibriiim  (adiabatic  mixing)  .  Note 
that  because  the  mixture  density  as  a  function  of  concentration  is  approxi¬ 
mated  by  adiabatic  mixing,  heat  transfer  is  not  included  in  the  DEGADIS 
simulations  of  the  DT  releases  herein. 

Table  10  shows  a  comparison  between  the  reported  temperature  associated 
with  the  maximum  reported  concentration  and  the  adiabatic  mixing  tempera¬ 
ture  associated  with  the  maximum  reported  concentration  for  the  DT  tests. 
The  assximption  of  no  heat  transfer  (or  adiabatic  mixing  of  the  released 
aerosol  with  the  ambient  air)  reproduces  the  reported  temperature  remark 
ably  well.  Furthermore,  Table  10  shows  the  difference  in  the  mixture 
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density  associated  with  the  difference  in  temperature  to  be  insignificant. 
Therefore,  the  assumption  of  insignificant  heat  transfer  is  justifiable  for 
the  DT  series.  The  validity  of  this  assumption  for  the  release  of  a  given 
material  is  generally  based  on  the  competing  effects  of  warming  the  con¬ 
taminant  mixture  by  mixing  with  (warm)  air  compared  with  warming  by  heat 
transfer:  both  of  these  effects  are  related  to  the  windspeed.  For  high 
windspeeds,  the  rate  of  mixing  is  increased,  thereby  adding  more  energy  to 
the  mixture.  Also  for  high  windspeeds,  the  mean  cloud  advection  speed  is 
higher,  implying  that  the  travel  time  to  a  given  distance  is  decreased,  and 
shorter  travel  times  provide  less  time  for  heat  transfer  to  occur.  It  is 
straightforward  to  show  that  the  mean  advection  time  to  a  given  distance 
should  be  proportional  to  (Tqq^  -  T^)  where  Tq^^  is  the  observed  tempera¬ 
ture  and  T^  is  the  adiabatic  mixing  temperature  for  the  same  concentration 
level.  For  lower  windspeeds  such  as  DT2  and  DT4,  the  heat  transfer  pro¬ 
cesses  were  more  important  than  for  DTI  and  DT3  which  had  higher  wind- 
speeds.  Indeed,  for  calm  conditions,  heat  transfer  processes  would  be  more 
important  than  the  heat  transfer  processes  observed  in  the  DT  series. 


TABLE  10.  COMPARISON  OF  THE  TEMPERATURE  ASSOCIATED  WITH  THE  OBSERVED 
MAXIMUM  CONCENTRATION  AND  THE  ADIABATIC  MIXING  TEMPERATURE 
ASSOCIATED  WITH  THE  MAXIMUM  CONCENTRATION  AT  800  METERS 


Measured 
Temperature 
at  the 

Maximum  Reported 
Concentration  ("C) 

Adiabatic  Mixing 
Temperature 
at  the 

Maximum  Reported 
Concentration  ("C) 

Density  Error 
in  Assuming  No 
Heat  Transfer  (%) 

DTI 

23.6 

20.9 

0.9 

DT2 

22.5 

17.2 

1.8 

DT3 

24.9 

22.4 

0.8 

DT4 

23.0 

17.7 

1.8 

The  exposure  limits  for  ammonia  (and  other  toxic  materials)  are 
generally  based  on  exposure  to  a  given  concentration  over  a  specified  time 
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period:  the  short-term  exposure  limit  (STEL,  based  on  a  15-m::.nute  exposure) 
for  ammonia  established  by  the  1979  American  Conference  of  Governmental 
Industrial  Hygienists  (ACGIH)  Is  35  parts  per  million  (Reference  20) . 
Therefore,  It  Is  necessary  to  be  able  to  predict  not  only  the  maximum 
concentration,  but  also  the  average  concentration  over  a  period  of  time 
(for  example,  15  minutes).  A  discussion  of  the  Implications  to  DEGADIS  of 
the  averaging  time  Is  included  in  Appendix  D.  In  the  DT  series,  the  MSA 
nondispersive  IR  gas  sensor  response  time  was  about  2  seconds,  and  the 
concentration  data  were  smoothed  using  a  3-second  sliding  average  (Refe¬ 
rence  1).  An  averaging  time  of  3  seconds  was  used  in  the  DEGADIS  simula¬ 
tions  for  comparison  with  the  data. 

The  Pasquill-Hanna  Gaussian  model  results  for  the  DT  series  are  shown 
In  Table  11  and  compared  with  the  concentration  profile  parameters  at  800 
meters  as  presented  in  Table  7.  As  indicated,  the  predicted  values  of 
are  In  reasonable  agreement  with  the  estimated  values,  but  the  predicted 
values  of  are  much  larger  (by  an  average  factor  of  about  6)  than  the 
estimated  values.  In  Tables  12-15  and  Figures  3-6,  the  model  results  are 
compared  with  maximum  reported  concentrations  (for  all  time) .  (All  model 
simulations  were  made  for  Che  ammonia  release  race  (Table  1)  and  the 
estimated  ammonia  mass  rate  at  800  meters  (Table  2).)  In  Tables  12-15,  the 
observed  temperatures  reported  correspond  In  time  with  the  maximum  reported 
concentrations.  The  values  of  Sy  and  are  reported  in  Tables  12-15;  the 
approximations  S  and  Sy  S  are  used  for  comparison.  As 

Indicated,  the  predicted  values  of  the  maximum  concentration  are  much 
smaller  (by  an  average  factor  of  about  6)  than  the  maximum  observed  values. 
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TABLE  11.  COMPARISON  OF  OBSERVED  AND  PREDICTED  GAUSSIAN 


PROFILE 

PARAMETERS 

AT  800  METERS 

Estimated 

Predicted* 

Estimated 

Predicted'* 

(m) 

a 

(a) 

(m) 

(m) 

DTI 

31.2 

31.2 

3.42 

21.0 

DT2 

65.7 

31.2 

3.60 

21.0 

DT3 

66.0 

31.2 

2.50 

21.0 

DT4 

86.7 

23.4 

2.80 

12.2 

a 

b 


From  Hanna  et  al.  (Reference  13)  and  corrected  for 
averaging  time 

From  Pasquill  and  Smith  (Reference  9) 


TABLE  12.  COMPARISON  OF  OBSERVED  AND  MODEL- PREDICTED  CONCENTRATION 
PROFILE  PARAMETERS  FOR  DTI 


at  3.5  km 

at  5.5  km 

Maximxim 
Concentra¬ 
tion 
(Vol.  %) 

Tempera¬ 

ture 

(°C) 

(m)  (m) 

Maximum 

Concentra¬ 

tion 

(Vol .  ppm) 

Maximum 

Concentra¬ 

tion 

(Vol .  ppm) 

Observed 

1.12 

23.6 

3.78  43.8 

670 

150 

Gaussian  Plume 

Model 

(66,5  kg/s) 

0.10 

-- 

21*  44^* 

75 

35 

(81  kg/s) 

0.12 

-- 

21*  44^* 

92 

43 

DEGADIS 

(66.5  kg/s) 

0.69 

24.0 

6.71  87.9 

850 

360 

(81  kg/s) 

0.83 

23.0 

6.06  90.5 

1040 

450 

*For  comparison,  S^  ■  a, 
^For  comparison,  Sy  3  Jla^ 
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TABLE  13.  COMPARISON  OF  OBSERVED  AND  MODEL- PREDICTED  CONCENTRATION 
PROFILE  PARAMETERS  FOR  DT2 


at  800  m 

at  1.4  km 

Maximum 
Concentra¬ 
tion 
(Vol.  %) 

Tempera¬ 

ture 

(“C) 

Sz 

(m) 

(m) 

Maximum 

Concentra¬ 

tion 

(Vol.  ppm) 

Observed 

1.86 

22.5 

4.80 

105 

4000 

Gaussian  Plume  Model 

(104.6  kg/s) 

0.21 

-- 

21* 

44b 

780 

(117  kg/s) 

0.24 

-- 

21* 

44b 

870 

DEGADIS 

(104.6  kg/s) 

1.59 

19.2 

3.53 

109 

5600 

(117  kg/s) 

1.81 

17.8 

3.36 

111 

6200 

®For  comparison,  S^  ” 
^For  comparison,  Sy  S 

TABLE  14.  COMPARISON  OF  OBSERVED  AND 
PROFILE  PARAMETERS  FOR  DT3 

MODEL- PREDICTED 

CONCENTRATION 

at  800  m 

at  2.8  km 

Maximum 
Concentra¬ 
tion 
(Vol.  %) 

Tempera¬ 

ture 

(”C) 

Sz 

(m) 

(m) 

Maximum 

Concentra¬ 

tion 

(Vol.  ppm) 

Observed 

1.62 

24.9 

2.66 

106 

1100 

Gaussian  Plume  Model 

(56.0  kg/s) 

0.08 

-- 

21* 

44b 

96 

(133  kg/s) 

0.20 

-- 

21* 

44b 

230 

DEGADIS 

(56.0  kg/s) 

0.61 

29.5 

7.30 

84.7 

1090 

(133  kg/s) 

1.37 

24.2 

4.82 

97.3 

2050 

®For  comparison,  S^  “  o- 
'^For  comparison,  Sy  • 
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TABLE  15.  COMPARISON  OF  OBSERVED  AND  MODEL- PREDICTED  CONCENTRATION 
PROFILE  PARAMETERS  FOR  DT4 


at  800  m 

at  2.8  km 

Maximum 
Concentra¬ 
tion 
(Vol.  %) 

Tempera¬ 

ture 

(“C) 

Sa 

(m) 

Sy 

(m) 

Maximum 

Concentra¬ 

tion 

(Vol .  ppm) 

Observed 

2.1 

23.0 

3.11 

143.5 

5100 

Gaussian  Plume  Model 

(72.1  kg/s) 

0.44 

-- 

12* 

33^ 

500 

(108  kg/s) 

0.65 

-- 

12* 

33^* 

750 

DEGADIS 

(72.1  kg/s) 

1.42 

22.7 

3.13 

97.5 

1900 

(108  kg/s) 

2.30 

17.0 

2.60 

102 

2600 

®For  comparison,  S^  ■  o_ 
^For  comparison,  Sy  2 


DEGADIS  predictions  (using  the  approximations  and  assumptions  described 
above)  for  the  DT  tests  are  shown  In  Tables  12-15  and  Figures  3-6.  The 
DEGADIS -predicted  maxlm\im  concentration  and  vertical  dispersion  parameter 
S2  are  generally  consistent  with  the  observed  values  at  800  meters.  It 
should  be  noted  that  the  maximum  reported  concentrations  for  the  sensor 
locations  past  1.4  kilometers  may  be  significantly  lover  than  the  maximum 
concentration  which  occurred  during  the  tests  due  to  the  wide  spacing  of 
these  sensors.  (Furthermore,  the  wind  field  may  not  have  been  constant  to 
these  long  distances . ) 

The  DEGADIS-predlcted  plume  width  (see  Figure  2)  Is  larger  than  the 
observed  width  due  to  the  presence  of  the  horizontally  homogeneous  central 
section  even  though  the  predicted  and  estimated  values  of  Sy  are  In  good 
agreement.  For  Test  4,  the  observed  width  to  the  0.2  percent  concentration 
level  was  440  meters  (Reference  21),  while  the  DEGADIS-predlcted  width  was 
800  meters  and  the  Gaussian  plume  model  prediction  was  60  meters  to  the 
same  concentration  level. 
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Maximum  Concentration  (%) 


Maximum  Concentration  (%) 


□  Experimental  Data 


Figure  A 


Distance  (m) 


Maximum  Observed  Concentration  and  Maximum  Predicted 
Concentrations  using  DEGADIS  and  the  Pasquill-Hanna 
Gaussian  Plume  Model  for  DT2.  (The  upper  line  for 
each  model  was  made  using  the  ammonia  mass  release  rate 
(Table  1),  and  the  lower  line  for  each  model  was  made 
using  the  estimated  ammonia  mass  rate  at  800  meters 
(Table  2). ) 


□  Experimental  Data 


1C?  1C?  1C? 

Distance  (m) 


Figure  5.  Maximum  Observed  Concentration  and  Maximum  Predicted 
Concentrations  using  DEGADIS  and  the  Pasquill-Hanna 
Gaussian  Plume  Model  for  DT3.  (The  upper  line  for 
each  model  was  made  using  the  ammonia  mass  release  rate 
(Table  1),  and  the  lower  line  for  each  model  was  made 
using  the  estimated  ammonia  mass  rate  at  800  meters 
(Table  2).) 


Maximum  Concentration  (%) 


Figure  6 


10^  ■  < - 1  I  I  i  HM - 1-  |.  I 

1^  1^  # 


Distance  (m) 


.  Maximum  Observed  Concentration  and  Maximum  Predicted 
Concentrations  using  DEGADIS  and  the  Pasquill-Hanna 
Gaussian  Plume  Model  for  DTA.  (The  upper  line  for 
each  model  was  made  using  the  ammonia  mass  release  rate 
(Table  1),  and  the  lower  line  for  each  model  was  made 
using  the  estimated  ammonia  mass  rate  at  800  meters 
(Table  2).) 


SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

Field-scale  releases  of  pressurized  liquid  anhydrous  ammonia  (the 
Desert  Tortoise  (DT)  series)  were  performed  by  Lawrence  Livermore  National 
Laboratories  (LLNL)  in  1983  for  the  U.S.  Coast  Guard,  the  Fertilizer 
Institute,  and  Environment  Canada.  Release  rates  for  the  four  experiments 
ranged  from  80  to  130  kilograms  per  second.  Ammonia  concentration  measure¬ 
ments  were  taken  at  sensor  arrays  located  100  and  800  meters  downwind  of 
the  release;  in  addition  to  these  arrays,  portable  gas  sensors  were 
deployed  from  1.4  to  5.5  kilometers  downwind  of  the  release.  The  data  from 
these  experiments  were  analyzed  to  determine  the  concentration  profile 
parameters  at  the  800 -meter  array.  These  concentration  profile  parameters 
were  compared  with  DEGADIS  and  the  Fasquill -Hanna  Gaussian  plume  model. 

The  Fasquill -Hanna  Gaussian  plume  model -predicted  maximum  concentrations 
were  significantly  lower  (by  an  average  factor  of  about  6)  than  the  maxlm\im 
reported  concentrations;  this  disagreement  is  due  to  the  overprediction  of 
the  vertical  dispersion  present  during  the  tests.  On  the  other  hand, 
DEGADIS -predicted  maximum  concentrations  were  in  reasonable  agreement  with 
the  maximum  reported  concentrations.  In  contrast  to  the  Gaussian  pl\une 
model,  DEGADIS  accounts  for  the  reduced  vertical  dispersion  due  to  the 
stable  density  stratifications  present  in  these  releases. 

A  quantitative  method  of  assessing  the  relative  important  of  Jet  and 
nonjet  dispersion  processes  was  developed.  Vhen  applied  to  the  DT 
releases,  the  jet  dispersion  processes  clearly  dominated  only  one  release 
(DT3)  according  to  the  developed  criteria. 

The  effects  of  different  averaging  times  on  the  DEGADIS  predictions 
were  discussed.  For  the  DT  series,  a  3-second  sliding  average  was  used  by 
LLNL  to  remove  noise  in  the  raw  data;  an  averaging  time  of  3  seconds  was 
used  in  the  DEGADIS  simulations. 

Based  on  comparison  of  reported  temperature  corresponding  to  the 
maximum  reported  concentration  and  the  adiabatic  mixing  temperature 
corresponding  to  the  same  maximum  reported  concentration,  heat  transfer  was 
found  to  be  insignificant  in  the  DT  series. 
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A  computer  code,  TSAUMA,  ves  used  to  determine  the  ammonla/alr/water 
mixture  density  as  a  fimctlon  of  ammonia  mole  fraction  for  Input  to 
DEGADIS.  Based  on  comparisons  with  TRAUMA  and  DEGADIS  computations,  It  was 
fotmd  th’it  assuming  Ideal  solution  liquid  phase  behavior  adequately 
describe  the  ammonla/alr/vater  mixture  density  for  the  DT  release  condi¬ 
tions.  It  was  also  fo\md  that  assuming  a  pseudobinary  of  alr/water  and 
ammonia  does  not  adequately  describe  the  ammonla/alr/water  mixture  density 
based  on  comparisons  with  DEGADIS. 

Finally,  a  linear  relationship  between  ammonia  concentration  and 
mixture  density  was  used  In  DEGADIS;  the  endpoints  were  represented  by  the 
released  ammonia  aerosol  and  ambient  air.  Upon  comparison,  the  DEGADIS - 
predicted  distances  using  this  simplified  scheme  were  within  about  +20 
percent  of  the  DEGADIS-predlcted  distance  to  a  given  concentration  level 
using  the  density  specification  from  TRAUMA  (Figure  1) .  This  last  simpli¬ 
fication  would  be  expected  to  be  less  valid  as  heat  transfer  from  the 
ground  to  the  cloud  becomes  more  Important.  Where  appropriate,  this  last 
simplification  may  be  useful  as  a  screening  technique  for  preliminary 
hazard  assessment. 
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APPENDIX  A 

DERIVATION  OF  THE  RELEASE  RICHARDSON  NUMBER  CRITERIA 


The  purpose  of  this  appendix  is  to  show  the  method  used  to  determine 
whether  a  grotmd-level  jet  dominates  a  release  or  whether  the  other  non¬ 
jet  phases  of  a  ground-level  release  (Including  the  negative  buoyancy- 
dominated  dispersion  phase,  the  stably  stratified  shear  flow  phase,  and  the 
passive  dispersion  phase)  are  the  only  Important  dispersion  phases. 

From  Havens  and  Spicer  (Reference  A-1)  and  Spicer  and  Havens  (Reference 
A-2).  the  criteria  for  determining  which  phase  of  the  dispersion  process 
dominates  a  particular  ground-level  release  was  based  on  water  txinnel 
experiments  reported  by  Brltter  (Reference  A-3).  In  that  set  of  experi¬ 
ments,  a  salt/water  solution  was  released  at  floor  level,  and  the  lateral 
and  upwind  extent  of  the  plume  was  recorded  as  a  function  of  the  buoyancy 
length  scale  used  by  Brltter  (Lj  -  Qg(pg  -  p4)/(pgU^)  where  Q  Is  the  volu¬ 
metric  emission  rate  and  u  Is  the  ambient  velocity).  Analysis  of  Brltter's 
releases  showed  that  the  release  was  passive  from  the  source  when  Lg/D  ^ 
0.005;  the  release  was  dominated  by  the  negative  buoyancy  dispersion  phase 
when  Lg/D  ?  0.1.  Based  on  these  observations,  the  following  criteria  of  a 
release  Richardson  number  were  obtained: 

If  Rlj.  5  30  negative  buoyancy-dominated  phase 

If  1  ^  Rl^  ^  30  stably  stratified  shear  flow  phase  (A-1) 

If  Rig  Z  1  passive  dispersion  phase 

where  Rig  -  g(Pg  -  These  values  were  obtained  using  the 

reported  ratio  (u/\i^)  «  16  for  the  water  flume.  Note  that  the  length  scale 
corresponding  to  the  depth  of  the  layer  Is  approximated  by  H  -  Q/uD  for 
these  nonjet  releases. 

With  the  criteria  of  domination  of  a  particular  nonjet  dispersion  phase 
established,  the  question  of  whether  a  particular  release  Is  dominated  by 
jet  effects  can  now  be  determined.  For  this  analysis,  a  release  will  be 
considered  dominated  by  jet  effects  associated  with  the  release  If  the  rate 
of  air  entrainment  due  to  the  jet  dominates  the  rate  of  air  entrainment  due 
to  the  dominant  nonjet  phase  of  the  release. 

When  the  release  Richardson  number  Is  less  than  1 ,  the  release  will  be 
dominated  by  the  passive  dispersion  phase  In  the  absence  of  jet  effects. 
Cude  (Reference  A-4)  and  Wheatley  (Reference  A-5)  report  that  for  a  passive 
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release,  jet  effects  will  dominate  a  particular  release  when  (V/u)  ?  1 
where  V  Is  the  velocity  of  the  released  material.  This  criterion  can  also 
be  obtained  by  examining  the  rate  of  air  entrainment  as  presented  by 
Vheatley  (Reference  A-6) 


dM 
_ a 

dx 


(¥] 


'j'’. 


(A-2) 


which  has  been  written  for  a  jet  released  near  ground  level  without  drag  on 
Che  bottom  surface.  The  DEGAOIS  model's  vertical  rate  of  air  entrainment 
per  unit  width  for  the  stably  stratified  shear  flow  and  passive  dispersion 
phases  Is  given  by 


dx 


^<R1*) 


(A-3) 


where  ^(Rl«)  -  0.88  +  0.099  Rl«.  Using  an  effective  width  of  2)rR,  a 
typical  value  ofa-0.2,  5l“2.1,  k-0.35,  and  (u/u*)  -  30  (typical  of 
atmospheric  boundary  layers) ,  the  above  equations  can  be  combined  to  show 
that  jet  effects  dominate  the  rate  of  dispersion  when 

-  >  16/(19  +  R1  )  (A-4) 

u  c 

This  criterion  has  the  characteristic  that  as  the  density  of  the  released 
fluid  Increases,  the  release  Is  more  readily  dominated  by  jet  effects 
(since  the  rate  of  air  entrainment  would  decrease  In  the  absence  of  the  jet 
effects).  Note  that  when  Rl^,  5!  1  (l.e.  a  passive  release),  this  criterion 
shows  that  the  Jet  effects  dominate  the  passive  dispersion  regime  when 
(V/u)  5  0.8,  which  Is  consistent  with  the  previously  reported  criterion. 

If  the  negative  buoyancy- dominated  dispersion  phase  Is  considered,  the 
ambient  flow  Is  no  longer  Important  In  determining  the  rate  of  non jet  air 
entrainment.  For  this  case,  the  relative  Importance  of  jet  effects  were 
evaluated  using  the  criterion  from  Britter's  data;  the  release  velocity  V 
was  used  in  place  of  the  ambient  velocity  In  the  buoyancy  length  scale.  It 
can  be  shown  that  jet  effects  dominate  the  negative  buoyancy -dominated 
phase  when 
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2 

>  10  Rig  (A- 5) 

This  criterion  has  the  characteristic  that  as  the  density  of  the  released 
fluid  increases,  the  release  is  less  readily  dominated  by  jet  effects 
(since  the  rate  of  air  entrainment  would  Increase  in  the  absence  of  the  jet 
effects  due  to  the  entrainment  associated  with  the  gravity  front) . 

The  following  procedure  is  then  suggested  for  determining  which  disper¬ 
sion  phase  is  dominant  from  the  start  of  a  release: 

(1)  Calculate  Rig  -  g(pg  -  Pg)H/(p^u£) . 

(2)  Determine  the  dominant  nonjet  dispersion  phase  in  the 
absence  of  a  jet  using  Equation  (A-1). 

(3)  Determine  if  ground-level  jet  effects  dominate  the  dominant 
nonjet  phase  determined  from  (2)  using  the  relationships 
summarized  in  Table  A-1. 

TABLE  A-1.  CRITERIA  FOR  DETERMINING  WHETHER  JET  EFFECTS  DOMINATE 
A  GROUND -LEVEL  RELEASE 


Ground- level  let  effects  dominate; 

Negative  buoyancy- dominated  phase  when 

Stably  stratified  shear  flow  phase  when 
Passive  dispersion  phase  when 


—  >  10  Ri 

I"*  j 

V/u  >  16/(19  +  Ri^) 
V/u  >0.8 
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APPENDIX  B 
DESCRIPTION  OF  THE 
DEGADIS  DENSE  GAS  DISPERSION  MODEL 


The  DEGADIS  (DEnse  GAs  DISperslon)  model  was  developed  from  research 
sponsored  by  the  U.S.  Coast  Guard  and  the  Gas  Research  Institute  (Reference 
B-1).  DEGADIS  Is  an  adaptation  of  the  Shell  HEGADAS  model  described  by 
Colenbrander  (Reference  B-2)  and  Colenbrander  and  Puttock  (Reference  B-3). 
DEGADIS  also  Incorporates  some  techniques  used  by  van  Ulden  (Reference 
B-4)  . 

If  the  primary  source  (gas)  release  rate  exceeds  the  maximTun  atmo¬ 
spheric  takeup  rate,  a  denser-than-alr  gas  blanket  Is  formed  over  the 
primary  source.  This  near-fleld,  buoyancy-dominated  regime  Is  modeled 
using  a  lumped  parameter  model  of  a  denser-than-alr  gas  "secondary  source" 
cloud  which  Incorporates  air  entrainment  at  the  gravity- spreading  front 
using  a  frontal  entrainment  velocity.  If  the  primary  source  release  rate 
does  not  exceed  the  maximum  atmospheric  takeup  rate,  the  released  gas  is 
taken  up  directly  by  the  atmosphere  and  dispersed  downwind.  For  either 
source  condition,  the  downwind  dispersion  phase  of  the  calculation  assumes 
a  power  law  concentration  distribution  In  the  vertical  direction  and  a 
modified  Gaussian  profile  in  the  horizontal  direction  with  a  power  law 
specification  for  the  wind  profile  (Figure  B-1).  The  source  model  repre¬ 
sents  a  spatially  averaged  concentration  of  gas  present  over  the  primary 
source,  while  the  downwind  dispersion  phase  of  the  calculation  models  an 
ensemble  average  of  the  concentration  downwind  of  the  source . 

A.  DENSER-THAN-AIR  GAS  SOURCE  CLOUD  FORMATION 

A  lumped  parameter  model  of  the  formation  of  the  denser-than-alr  gas 
source  cloud  or  blanket,  which  may  be  formed  from  a  primary  source  such  as 
an  evaporating  liquid  pool  or  otherwise  specified  ground-level  emission 
source,  or  by  an  initially  specified  gas  volume  of  prescribed  dimensions 
for  an  Instantaneous  release,  is  Illustrated  in  Figure  B-1.  The  gas 
blanket  is  represented  as  a  cylindrical  gas  volume  which  spreads  laterally 
as  a  density -driven  flow  with  entrainment  from  the  top  of  the  source 
blanket  by  wind  shear  and  air  entrainment  into  the  advancing  front  edge. 

The  source  blanket  will  continue  to  grow  over  the  primary  source  until 
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Schematic  Diagram  of  DEGADIS  Dense  Gas  Dispersion  Model  (Concluded) 


the  atmospheric  takeup  rate  from  the  top  Is  matched  by  the  air  entrainment 
rate  from  the  side  and,  If  applicable,  by  the  rate  of  gas  addition  from 
under  the  blanket.  Of  course,  the  blanket  is  not  formed  if  the  atmospheric 
takeup  rate  is  greater  than  the  evolution  rate  of  the  primary  source .  For 
application  of  the  downwind  calculation  procedure,  the  blanket  is  modeled 
as  being  stationary  over  the  center  of  the  source  (x  -  0) . 

1.  Secondary  Source  Blanket  Extent  for  Ground-Level  Releases 

If  a  denser -than -air  gas  blanket  is  present,  the  (downwind)  emis¬ 
sion  rate  from  the  blanket  is  equal  to  the  maximum  atmospheric  takeup  rate. 

O 

That  is,  for  E(t)/jrR^(t)  >  *  source  blanket  is  formed  over  the 

primary  source.  The  blanket  frontal  (spreading)  velocity  is  modeled  as 


f  E 


(B-1) 


where  p  is  the  average  density  of  the  source  blanket.  This  gravity 
intrusion  relationship  is  applicable  only  for  p  >  the  value  of  C£  used 
is  1.15  based  on  laboratory  measurements  of  cloud  spreading  velocity 
(Reference  B-1). 

The  blanket  radius  R  as  a  function  of  time  is  determined  by 
integrating  dR/dt  -  U£.  When  the  total  mass  of  the  cloud  is  decreasing 
with  time,  the  radius  is  assumed  to  decrease  according  to  (dR/dt)/R 
-  (dH/dt)/H  for  ground-level  sources.  The  radius  of  the  blanket  is 
constrained  to  be  greater  than  or  equal  to  the  radius  Rp  of  any  primary 
(liquid)  source  present. 


2.  Secondary  Source  Blanket  Extent  for  Instantaneous  Releases 
The  gravity  intrusion  relationship  (Eq\iation  (B-1))  will 
overpredict  initial  velocities  for  instantaneous,  abovegro\ind  releases  of  a 
denser-than-air  gas  since  no  initial  acceleration  phase  is  included.  In 
this  case,  the  following  procedure  adapted  from  van  Ulden  (Reference  B-4) 
is  recommended. 
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For  Instantaneous  gas  releases,  the  radially  S3nnmetric  cloud  is 
considered  to  be  composed  of  a  tall  section  with  height  and  radius 
and  a  head  section  with  height  (Figure  B-2) .  A  momentum  balance  is  used 
to  accoTxnt  for  the  acceleration  of  the  cloud  from  rest;  the  effect  of 
ambient  (wind)  momentum  Is  Ignored.  Although  the  following  equations  are 
derived  asstaming  the  prlmairy  source  emission  rate  Is  zero,  the  resulting 
equations  are  assumed  to  model  the  secondary  source  cloud  development  when 
the  primary  source  rate  Is  nonzero.  When  the  frontal  velocity  from  the 
momentum  balance  Is  the  same  as  Equation  (B-1),  the  momentu.m  balance  is  no 
longer  applied  and  the  frontal  velocity  is  given  by  Equation  (B-1) . 
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Figure  B-2.  Schematic  Diagram  of  a  Radially  Spreading  Cloud. 


Three  mein  forces  act  on  the  cloud:  (1)  a  static  pressure  force 
(Fp) ,  (2)  a  dynamic  drag  force  <F^),  and  (3)  a  force  which  accounts  for  the 
acceleration  reaction  of  the  ambient  fluid,  represented  as  a  rate  of 
virtual  momentum  change  with  respect  to  time  (-dP^/dt).  Denoting  the 
momentum  of  the  head  and  tall  as  and  P^,  respectively,  the  momentum 
balance  Is 


or 


dt 


dt 


(Ph  >  P,) 


F 

P 


^  (P.  +  P^  +  P  )  -  F  +  Fj 
dt  h  t  v'  p  d 


(B-2) 


(B-3) 


The  terms  In  the  momentum  balance  are  evaltiated  differently  for 
early  times  before  a  gravity  current  head  has  developed  (Hj^  <  H,.)  and  for 
times  after  the  head  has  developed  but  the  cloud  Is  still  accelerating 
(Figure  B-2).  Because  the  gravity  current  head  develops  so  rapidly,  the 
model  equations  describing  the  times  after  the  gravity  current  head  forms 
(Hh  ^  Hj.)  are  derived  first.  The  model  equations  describing  earlier  times 
(Hjj  <  H^)  use  simplification  of  the  equations  for  >  H^. 


a.  Unsteady  Gravity  Current 

When  the  cloud  accelerates  to  the  point  that  Hjj  >  (Figures 
B-2,  B-3),  the  frontal  velocity  Is  determined  from  the  momentum  balance 
(Equation  (B-2))  as  follows. 

The  static  pressure  force,  obtained  by  Integrating  the  static 
pressure  over  the  boundary  of  the  current.  Is 

Fp  -  j  i  gdpH^  j  J  2,rRH^  J  -  .rgApRH^  (B-4) 

Neglecting  the  shear  stress  at  the  bottom,  the  dynamic  force  on  the  current 
Is  the  sum  of  the  drag  force  on  the  head  of  the  current  and  the  lift 
forcethat  arises  due  to  asymmetry  In  the  ambient  flow  around  the  head.  The 
drag  force  Is  represented  by 

^  -  r  '."f  ( 2'Vv“h )  -  <“-5) 
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Figure  B-3.  The  Unsteady  Gravity  Current  (Reference  B-4) . 


where  Is  an  effective  drag  coefficient  and  the  constant  a^  Is  an 
empirical  ratio  of  the  average  head  depth  to  (a^  -  • 

The  horizontal  acceleration  reaction  (-dP^/dt)  Is  approximated 
by  the  reaction  to  an  accelerating  elliptical  cylinder  with  an  aspect  ratio 
H/R  (Reference  B-5): 

(  ‘'l^a’^“*S  )  (B-6) 

and  the  vertical  acceleration  reaction  Is  represented  as 

where  kj^  and  k2  are  coefficients  of  order  one.  Using  a  single  constant, 
Equations  (B-6)  and  (B-7)  give 


dP 

V 

dt 
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(B-8) 


i- 

t- 


i 


I 


dP^  d(RH^u^) 

’  dt  "  dt 

Using  Equations  (B-4) ,  (B-5),  and  (B-8),  the 
B-2))  becomes 

<iP  2  2 

-T-  -  jrgApRH^  -  a  d  irp  RH.  u-  -  e  irp 
dt  *  ^  t  V  V  '^a  ^  f  V  '^a 

Following  van  Olden  (Reference  B-4,  Reference  B-6),  it  is 
assumed  that  the  potential  energy  decrease  due  to  slumping  of  the  cloud  is 
offset  by  the  production  of  kinetic  energy,  which  through  the  action  of 
shear,  is  partly  transformed  to  turbulent  kinetic  energy.  Fart  of  the 
turbulent  kinetic  energy  is  transformed  back  into  potential  energy  due  to 
entrainment  of  air  by  the  cloud.  This  "buoyant  destruction"  of  kinetic 
energy  is  assumed  to  be  proportional  to  the  rate  of  production  of  turbulent 
kinetic  energy,  and  following  Simpson  and  Britter  (Reference  B-7)  it  is 
assumed  that  the  turbulent  kinetic  energy  production  rate  scales  as 
irp^HRu| .  Then , 

I  gApH  “  -  eirp^HRu^  (B-10) 


momentum  balance  (Equation 


d(RH  u^) 
dt 


(B-9) 


which  can  be  written 

e(2irRH)u. 


dt 


gApH 

2 


€(2jrRH)u. 


Ri- 


(B-11) 


where  e  is  an  empirically  determined  coefficient, 
the  air  entrainment  rate. 


Noting  dV/dt  represents 


(B-12) 


where  represents  the  air  entrainment  mass  rate. 
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The  volume  Integral 


V  -  2w  r  h(r,t)rdr 

Jq 

where  h(r,t)  is  to  be  expressed  in  terms  of  and  H 
momentum  Integral 

P  -  2w  I  pu<r,t)h<r,t)rdr  -  +  P. 

Jq  t  h 


t’ 


and  the 


(B-13) 


(B-14) 


are  then  approximated  with  separate  analyses  of  the  head  and  tall  of  the 
current . 

In  the  tail  of  the  current,  the  shallow  water  equations  are 
assumed  applicable.  It  is  assumed  that  the  shape  of  the  current  is  quasi- 
statlonary  in  time,  and  the  layer-averaged  density  difference  is  assumed 
horizontally  uniform.  It  follows  that  the  volume  and  momentum  of  the  tail 
are  given  by 

’t  -  (  V  “h  )  /  O-is) 


^  - 1 '  ( I  "t  *  'S, )  r 

A  momentum  balance  for  the  head  region,  Figure  B-4,  assuming 
quasi-steady  state,  Indicates  that  the  static  and  dynamic  pressure  forces 
on  the  head  should  be  balanced  by  the  net  fl\ix  of  momenttim  due  to  flow  into 
and  out  of  the  head.  The  static  pressure  and  drag  are,  respectively 

Fp  -  (  i  (B-17) 

'C  -■'*»(  I  '."J  ]  [  ] 


Near  the  surface,  the  inward  flow  (u^  in  Figure  B-4)  carries  momentum  into 
the  head,  while  the  return  flow  <U3  in  Figure  B-4)  carries  momentum  out  of 


Figure  B-4.  The  Heed  of  a  Steady  Gravity  Current 
(Reference  B-4,  Reference  B-7). 


the  head.  Asstuning  uj  =»  u^.  =  1/2  and  u^  =»  ^v'^f'  momentum  flux 

Into  the  head  Is  approximately 

\  *  [  2'»h“h  ] 

Upon  rearranging,  the  momentum  balance  on  the  head  gives 

^  .  l./(  .A  -  ]  -  <i  (B-20) 

when  »  0.2  and  -  0.64;  Equation  (B-20)  then  specifies  the  head 
velocity  boundary  condition.  The  volume  of  the  head  is  determined  by 
assuming  that  the  head  length  scales  with  It  follows  that 

R  -  Rjj  “  Vl 

where  b.^  Is  an  empirical  constant,  and  the  volume  of  the  head  becomes 

If  the  layer -averaged  velocity  Is  assumed  to  increase  linearly  with  r,  it 
follows  that 
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"h-"f 


(B-23) 


R 


and 


2. 

3  R 


Along  with  the  definition  of  U£, 


dt 


—  u 


f  ’ 


(B-24) 


(B-25) 


Equations  (B-9).  (B-11).  (B-20).  (B-21) ,  (B-23).  and  (B-25)  are  solved  to 
determine  p,  Hjj,  V,  and  when 

The  constants  a^,  b^,  d^,  e^,  and  t  are  assigned  values  1.3, 
1.2,  0.64,  20.,  and  0.59,  respectively,  based  on  analysis  of  the  stlll-alr 
denser-than-alr  gas  release  experiments  of  Havens  and  Spicer  (Reference 
B-1). 


b.  Initial  Gravity  Current  Development 

To  model  the  Initial  cloud  shape,  the  tail  and  head  height  are 
considered  constant  with  respect  to  radius.  The  momentum  balance  on  the 
cloud  Is  then  given  by 


dt 


jrgAp 


+  a  b  HT 
V  V  n 


dP 


-  ’•vV."S.“£  -  dT 


(B-26) 


where  the  first  term  on  the  right-hand  side  represents  the  static  pressure 
force  on  the  head  and  the  second  term  represents  the  drag  force  on  the 
bottom  surface  of  the  cloud.  The  third  force  is  the  acceleration  reaction 
by  the  ambient  fluid,  represented  by  Equation  (B-8) . 

The  dimensions  of  the  head  are  again  given  by 

\  -  R  -  \\\  <B-27> 
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and 


u  2 


[gAp/P^] 


(B-28) 


When  the  height  of  the  tail  Is  assumed  uniform  with  respect  to  radius , 
it  follows  that 


“t  -  !  -  "Ik  *  *1.) 


2 

(xR^) 


(B-29) 


where  M  is  the  total  mass  of  the  cloud.  The  momentum  of  the  head  and 
tail  are  then 


2 

\  ■  3  **v - ;; - “f 


(B-30) 


2 


^  *  3 '  “r"  '"f 


(B-31) 


Equations  <B-26)  through  <B-31)  determine  the  momentum  of  the  blanket  as  a 
function  of  time,  and  thus  the  frontal  velocity  uj.  The  cloud  accelerates 
from  rest  because  »  0  initially. 


3 .  Material  and  Energy  Balances 

The  balance  on  the  total  mass  of  gas  in  the  source  blanket 
(M  -  jrR^Hp)  is 


dt  “  dt 


E<t)  +  M  +  M 

a  w,s 


f  1  (-1 


(B-32) 


where  E(t)  is  the  gas  evolution  rate  from  the  primary  (liquid)  source. 
For  spills  over  water,  the  water  entrainment  term  (M  )  is  included  in 
the  source  blanket  description  and  is  calculated  from  Equation  (B-46), 
and  the  (humid)  air  entrainment  rate  (Equation  (B-12))  is 


-  2»rRH(£U-)p^ 
a  fa 


gApH/(p^u^)  j 


(B-33) 


U8 


The  balance  on  the  mass  of  contaminant  in  the  source  blanket 


2 

(M  —  w  jtR  Hp)  is 
c  c 


dT  -  ] 


and  the  mass  balance  on  the  air  In  the  source  blanket  (M  -  v  ttR  Hp) 

a  a 


dM  . 
_ a  ^ 

dt  "  dt 


w  irR  Hp 
a 


1  M 

1-r^ 


*max 


w  (»rR‘) 
a 


(B-34) 

is 


(B-35) 


where  the  ambient  humidity  is  and  the  mass  fraction  of  contaminant  and 
air  are  w^.  »  M^/M  and  -  M^/M,  respectively. 
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The  energy  balance  on  the  source  blanket  (hjrR  Hp)  gives 


dt 


hjrR^Hp 


h_E(t)  +  h  M  +  h  M 
E  a  a  w  v,s 


♦max 


(♦R  )  +  Q. 


(B-36) 


where  hg  is  the  enthalpy  of  the  emitted  gas,  h^  is  the  enthalpy  of  the 
ambient  hrimid  air,  and  h^  is  the  enthalpy  of  any  water  vapor 
entrained  by  the  blanket  if  over  water.  There  are  three  alternate 
submodels  included  for  the  heat  transfer  (Q^)  from  the  surface  to  the 
c  lotid 

The  simplest  method  for  calculating  the  heat  transfer  between  the 
substrate  and  the  gas  cloud  is  to  specify  a  constant  heat  transfer 
coefficient  for  the  heat  transfer  relation 

Q,  -  q,  [  '  (  -  »p  ]  ]  -  -  "p  )  ] 

where  is  the  rate  of  heat  transfer  to  the  cloud,  is  the  heat  fltix, 
and  AT  is  the  temperature  difference.  For  the  calculation  of  heat  transfer 
over  the  source ,  the  temperature  difference  is  based  on  the  average  tem¬ 
perature  of  the  blanket. 


In  the  evaluation  of  the  Burro  and  Coyote  series  of  experiments, 
Koopman  et  al.  (Reference  B-8)  proposed  the  following  empirical  heat  trans¬ 
fer  coefficient  relationship  for  heat  transfer  between  a  cold  LNG  cloud  and 
the  ground 

ho  -  (B-38) 


where  the  value  of  Vy  was  estimated  to  be  0.0125  m/s.  This  constant  can  be 
varied  in  the  model . 

From  the  heat  transfer  coefficient  descriptions  for  heat  transfer 
from  a  flat  plate,  the  following  relationships  can  be  applied.  For  natural 
convection,  the  heat  transfer  coefficient  is  estimated  using  the  Nusselt 
(Nu) ,  Grashoff  (Gr) ,  and  Schmidt  (Sc)  ntimbers  (Reference  B-9)  from 

Nu  -  0.14  (Gr  Sc)^^^  (B-39) 


Or 


h  -  0.14 
n 


2^3 

gP  C  « 

- ^  AT 

T  Pr 


1/3 


(B-40) 


where  h^^  is  the  heat  transfer  coefficient  due  to  natural  convection  and  Pr 
is  the  Prandtl  number.  To  simplify  the  calculations,  the  parameter  group 


Pr 


-2 


( “p™  )M  ] 


1/3 


(B-41) 


is  estimated  to  be  60  in  mks  units.  The  actual  value  of  the  group  is 
47.25,  58.5,  and  73.4  for  air,  methane,  and  propane,  respectively. 
Equation  (B-40)  becomes 


h  -  18 
n 


(B-42) 


where  the  density  p,  molecular  weight  MW,  and  temperature  difference  AT  are 
based  on  the  average  composition  of  the  gas  blanket. 
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For  forced  convection,  the  Colburn  analogy  (Reference  B-10)  Is 
applied  to  a- flat  plate  using  the  Stanton  number  for  heat  transfer  St^^  and 
the  Frandtl  number  as 


2/3 

ct  Pr  '  «  —  —  _ 

H  2  - 

u 


(B-43) 


-2/3 


h^  -  (upCp)  Pr 


where  h^  is  the  heat  transfer  coefficient  due  to  forced  convection, 
velocity  is  evalviated  at  z  =  H/2  and  Pr  is  estimated  to  be  0.741, 


(B-44) 


If  the 


(B-45) 


If  H/2  <  Zp^,  then  the  velocity  is  evaluated  at  z  >  Zp. 

The  overall  heat  transfer  coefficient  is  then  the  maximum  of  the 
forced  and  natural  coefficients,  i.e.  hg  -  max(h£,h^).  The  heat  flux  and 
transfer  rate  are  then  estimated  by  Equation  (B-37). 

If  the  gas  blanket  is  formed  over  water,  water  will  be  transferred 
from  the  surface  to  the  cloud  by  a  partial  pressure  driving  force 
associated  with  the  temperature  difference  between  the  surface  and  the  gas 
blanket.  The  rate  of  mass  transfer  of  water  is 


r  ( Cs  -  ».,c  ]  [  ’  ( 


2  2 
R  -  R^ 
P 


(B-46) 


where  Fq  Is  the  overall  mass  transfer  coefficient.  The  driving  force  is 
the  difference  of  the  vapor  pressure  of  water  at  the  surface  temperature 

'ft 

p^  g  and  the  partial  pressure  of  water  in  the  cloud,  p^  .  (The  water 
partial  pressure  in  the  cloud  is  the  minim\im  of:  (1)  the  water  mole  frac¬ 
tion  times  the  ambient  pressure;  or  (2)  the  water  vapor  pressure  at  the 

'ft 

cloud  temperature  (p^  c^'^  natural  convection  coefficient  is  based  on 

the  heat  transfer  coefficient  and  the  analogy  between  the  Sherwood  number 
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(Sh)  and  the  Nusselt  number  (Nu)  suggested  by  Bird  et  al.  (Reference  B-11) 


Sh  -  Nu  -  0.14  (Gr  - 

If  the  Schmidt  number  Is  taken  as  0.6 

.9 

2.2  X  10  in  mks  units. 


) 


(B-47) 

is  estimated  to  be 


F  -  9.9  X  10 
n 


(B-48) 


For  forced  convection.  Treybal  (Reference  B-10)  suggests  that  the  Stanton 
number  for  mass  transfer  Stjj  and  the  Stanton  number  for  heat  transfer  Sty 
are  related  by 


Or. 


r  Pr 

[  S^  J  -  ^-^5  Stjj 

20.7  h 

F  »  - - 

f  MW  C 

P 


(B-49) 

(B-50) 


The  overall  mass  transfer  coefficient  Fq  is  calculated  as  the  larger  of  the 
natural  and  forced  convection  coefficients. 

For  the  case  when  the  primary  (liquid)  source  emission  rate  E(t)  is 
larger  than  the  atmospheric  takeup  rate  ,  Equations  (B-32), 

(B-34),  (B-35),  and  (B-36)  are  integrated  for  the  mass,  concentration,  and 
enthalpy  of  the  gas  blanket  along  with  an  appropriate  equation  of  state 
(i.e.  relationship  between  enthalpy  and  temperature  and  between  temperature 
and  density). 

For  the  case  when  the  emission  rate  is  not  sufficient  to  form  a  gas 
blanket,  the  flux  of  contaminant  is  not  determined  by  the  maximum  atmo¬ 
spheric  takeup  rate.  Consider  the  boundary  layer  formed  by  the  emission  of 
gas  into  the  atmosphere  above  the  primary  source.  If  the  source  is  modeled 
to  have  a  uniform  width  2b  and  entrain  no  air  along  the  sides  of  the  layer, 
the  balance  on  the  total  material  (pj^Uj^Hj^)  in  a  differential  slice  of  the 
layer  is 


r  „  1 

f  Q* 

PtU-H. 

=•  p  w  + 
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(B-51) 
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where  w-  Is  the  vertical  rate  of  of  air  entrainment  into  the  layer  given  by 
Equation  (B-i3),  t**®  average  density  of  the  slice,  and  Q*/Wj,  is  the 

total  flux  of  gas  from  the  primary  (liquid)  source.  The  balance  on  the 
mass  flow  rate  of  contaminant  '  *up^ 

S.l'W  -  -  'up) 

With  an  equation  of  state  to  relate  Cg  ^  ^L*  Eqtiation  (B-51)  is  inte¬ 
grated  from  the  upwind  edge  of  the  source  (z  «  ^p)  downwind  edge 

(X  -  L  +  x^p) . 

In  order  to  generate  the  initial  conditions  for  the  downwind  dis¬ 
persion  calculations,  the  maximum  concentration  Cg  and  the  vertical  dis¬ 
persion  parameter  S2  are  needed.  Since  Equations  (B-Sl)  and  (B-52)  are 
written  for  a  vertically  averaged  layer,  consider  the  vertical  average  of 
the  power  law  distribution.  The  height  of  the  layer  is  the  height  to 
some  concentration  level,  say  10  percent  of  the  maximum.  Although  strictly 
a  function  of  a,  this  value  is  modeled  by 

\  -  ‘l  “eff 


where  Hgj-p  is  the  effective  height  defined  by  Equation  (B-79)  and  6^ 
2.15.  The  vertically  averaged  concentration  Cg  ^  can  be  defined  by 


A  -  r 


cdz 


(B-54) 


And  similarly,  the  effective  transport  velocity  u^^  is  defined  by 


-  r 


cu  dz 
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(B-55) 


With  Equation  (B-53)  and  defining  relations  for  Hgpp  and  Ugpp 
(Equations  (B-79)  and  (B-93),  respectively),  it  follows  that 
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and 
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(B-58) 


w 
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where  w^  Is  given  by  Equation  (B-83) . 


4.  Maximum  Atmospheric  Takeup  Rate 

The  maximum  atmospheric  takeup  rate  will  be  the  largest  takeup  rate 
which  satisfies  Equations  (B-51)  and  (B-52).  As  well,  the  maximum  concen¬ 
tration  of  contaminant  in  the  power  law  profile  at  the  downwind  edge  of  the 
source  will  be  the  source  contaminant  concentration  (Cj.)^.  If  Equations 
(B-51)  and  (B-52)  are  combined  along  with  the  assumption  of  adiabatic 
mixing  of  Ideal  gases  with  the  same  constant  molal  heat 


capacity  (l.e. 


7  -  constant) ,  the  maximum  takeup  flux  is 


modeled  by 
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*max 


c  s 


ku^d  +  a) 
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(B-59) 


where 


:  ”  L  Iq  ^(Ri^) 

where  <^(R1  )  is  given  in  Equation  (B-76)  for  p  >  p  . 

"  a 

An  upper  bound  of  the  atmospheric  takeup  flux  can  be  characterized 
by  the  condition  where  the  source  begins  to  spread  as  a  gravity  intrusion 
against  the  approach  flow.  In  water  flume  experiments,  Britter  (Reference 
B-12)  measured  the  upstream  and  lateral  extent  of  a  steady- state  plume  from 
a  circular  source  as  a  function  of  Rl*.  A  significant  upstream  spread  was 
obtained  for  Rl^  >  32,  and  lateral  spreading  at  the  center  of  the  source 
was  Insignificant  for  Rl*  <  8.  The  presence  of  any  significant  lateral 
spreading  represents  a  lower  boxind  on  the  conditions  of  the  maximum  takeup 
flux. 

The  integral  of  Equation  (B-60)  is  calculated  using  a  local 
Richardson  number  of 
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where 
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(B-61) 
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and  <f>^  is  3.1  (corresponding  to  Rt^  -  20(8  <  Ri^  <  32)).  Using  this 
Rl^(x) ,  Eqtiatlon  (B-60)  Is 


1  1  dx 

I  “  L  Jo  1.04 

^  0.88  +  0.099 


To  simplify  the  numerical  problem,  the  Integral  Is  approximated  as 
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which  then  specifies  the  maximum  atmospheric  takeup  fl\ix. 


5.  Transient  Denser-than-Alr  Gas  Release  Simulation 

If  a  steady- state  spill  Is  being  simulated,  the  transient  source 

calculation  Is  carried  out  until  the  source  characteristics  are  no  longer 

varying  significantly  with  time.  The  maximum  centerline  concentration  c^, 

the  horizontal  and  vertical  dispersion  parameters  S„  and  S_,  the  half  width 

y  * 

b,  and  If  necessary,  the  enthalpy  h  are  used  as  Initial  conditions  for  the 
downwind  calculation  specified  In  a  transient  spill. 

If  a  transient  spill  Is  being  simulated,  the  spill  Is  modeled  as  a 
series  of  pseudo-steady-state  releases.  Consider  a  series  of  observers 
traveling  with  the  wind  over  the  transient  gas  source  described  above;  each 
observer  originates  from  the  point  which  corresponds  with  the  maximum 
upwind  extent  of  the  gas  blanket  (x  -  -R^ax^ •  desired  observer 
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velocity  Is  the  average  transport  velocity  of  the  gas  U£pp  from  Equation 
(B-93);  however,  the  value  of  u^pp  will  differ  from  observer  to  observer 
with  the  consequence  that  some  observers  may  be  overtaken  by  others .  For  a 
neutrally  buoyant  cloud,  Ugpp  becomes  a  ftmction  of  downwind  distance  alone 
which  circumvents  this  problem.  With  this  functionality,  Colenbrander 
(Reference  B-2)  models  the  observer  velocity  as 


(B-64) 


2 

where  S  is  the  value  of  S  when  the  averaged  source  rate  (jtR  Q  )  is 
^Om  ^0 

a  maximum  and  the  subscript  i  denotes  observer  i.  Noting  that 

Uj^(x)  -  dx^/dt,  observer  position  and  velocity  as  functions  of  time  are 

determined. 

A  pseudo-steady-state  approximation  of  the  transient  source  is 
obtained  as  each  observer  passes  over  the  source.  If  t  and  t, 

upj^  dnj^ 

denote  the  times  when  observer  1  encounters  the  upwind  and  downwind 
edges  of  the  source,  respectively,  then  the  source  fetch  seen  by  observer 
1  is; 


L,  -  X  -  X . 
i  upj^  dnj^ 

« 

The  width  of  the  source  2Bj^(t)  is  defined  by 
B'^t)  -  R^(t)  -  xj(t) 

Then  the  gas  source  area  seen  by  observer  i  is 

r^'dn 

2Lj^b  -  2 

uPi 


(B-65) 
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where  2b ^  is  the  average  width. 

The  takeup  rate  of  contaminant  2(Q^Lb)^  is  calculated  as 
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The  total  mass  flux  rate  from  the  source  is 

f 

t  c 


(B-69) 


Vlth  these  equations,  the  average  composition  of  the  layer  can 

be  determined  at  each  x  -  x  over  Che  source.  With  the  enthalpy 

up 

of  the  layer  given  by 
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(due  to  the  choice  of  the  reference  temperature  as  the  ambient  temperature) 
and  with  a  suitable  eq\iation  of  state  relating  enthalpy,  temperature,  and 
density,  Che  source  can  be  averaged  for  each  observer.  After  the  average 
composition  of  the  layer  is  determined  at  the  downwind  edge,  an  adiabatic 
mixing  calculation  is  performed  between  this  gas  and  the  ambient  air.  This 
calculation  represents  the  function  between  density  and  concentration  for 
the  remainder  of  the  calculation  if  the  calculation  is  adiabatic;  it  repre¬ 
sents  the  adiabatic  mixing  condition  if  heat  transfer  is  included  in  the 
downwind  calculation. 

For  each  of  several  observers  released  successively  from  x  -  ~^ax’ 
the  observed  dimensions  L  and  b,  the  downwind  edge  of  the  source  x^^,  the 
average  vertical  dispersion  coefficient  S^,  the  average  takeup  flux  Q«,  the 
centerline  concentration  c^.,  and  if  applicable,  the  average  enthalpy  h^  can 
be  determined  for  each  observer.  With  these  input  values,  a  steady- state 
calculation  is  made  for  each  observer.  The  distribution  parameters  for  any 
specified  time  t^  are  determined  by  locating  the  position  of  the  series  of 
observers  at  time  t^,  l.e.  x^Ct^).  The  corresponding  concentration  distri¬ 
bution  is  then  computed  from  the  assumed  profiles. 


B.  STEADY- STATE  DOWNWIND  DISPERSION 

The  model  treats  dispersion  of  gas  entrained  Into  the  vlnd  field  from 
an  idealised,  rectangularly  shaped  source  of  width  2b  and  length  L.  The 
circular  source  cloud  Is  represented  as  an  equivalent  area  rectangle 
(jtR^  -  2bL)  with  equivalent  fetch  (L  -  2R) .  Similarity  forms  for  the  con¬ 
centration  profiles  are  assumed  which  represent  the  plume  as  being  composed 
of  a  horizontally  homogeneous  section  with  Gaussian  concentration  profile 
edges  as  follows : 


c(x,y,z)  - 


S  (x) 

y 


S^(x) 


for  |y|  >  b 


-  c^<x)  exp 


S^(x) 


A  power  law  wind  velocity  profile  is  ssstimed 


u  “  u-  1  — 
X  O  s, 


for  |y|  <  b 
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(B-72) 


where  the  value  of  a  is  determined  by  a  weighted  least-squares  fit 
of  the  logarithmic  profile 


u  -  ; —  In 

X  k 


z  +  z. 


(B-73) 


Functional  forms  for  ^  and  typical  values  of  a  are  given  in  Table  B-1  for 
different  Fasquill  stability  categories.  With  these  profiles,  the 
parameters  of  Eqtiation  (B-71)  are  constrained  by  ordinary  differential 
equations. 


1.  Vertical  Dispersion 

The  vertical  dispersion  parameter  is  determined  by  requiring 
that  it  satisfy  the  diffusion  equation 


Curve  fit.  of  data  from  Pasquill  (Reference  17) 


(B-74) 


3c  _  ^  3c 
3x  3z  z  3z 


with  the  vertical  turbulent  dlffuslvlty  given  by 


*-z  ”  ^(Rl.) 


(B-75) 


The  function  ^(Rl^)  Is  a  curve  fit  of  laboratory- scale  data  for 
vertical  mixing  In  stably  denslty-stratlf led  fluid  flows  reported  by  Kantha 
et  al.  (Reference  B-2),  Lofqulst  (Reference  B-14),  and  McQuald  (Reference 
B-15)  for  Rl*  >  0.  For  Rl*  <  0,  the  function  ^(Rl^)  Is  taken  from  Colen- 
brander  and  Puttock  (Reference  B-3)  and  has  been  modified  so  the  passive 
limit  of  the  two  functions  agree  as  follows: 


.^(Rl.)  -  0.88  +  0.099  Rr 
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The  friction  velocity  Is  calculated  using  Eqxiatlon  (B-73)  from  a 
known  velocity  uq  at  a  specific  height  zq.  Combining  the  assxmed  simi¬ 
larity  forms  for  concentration  and  velocity,  Eqtiatlons  (B-71),  (B-72), 
(B-74).  and  (B-75)  give 
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where  the  Richardson  number  Rl  Is  computed  as 
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and  the  effective  cloud  depth  Is  defined  as 
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Equation  (B-77)  can  be  viewed  as  a  volumetric  balance  on  a  differ¬ 
ential  slice  of  material  downwind  of  the  source.  For  a  mass  balance  over 
the  same  slice, 

5  (  'lVl  ]  - 

which  Is  the  same  result  as  Equation  (B-51)  without  the  source  term. 

With  Equations  (B-57)  and  (B-58),  this  becomes 
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Using  the  assumption  of  adiabatic  mixing  of  Ideal  gases  with  the 

same  constant  molal  heat  capacity  (l.e.  - ^  -  constant)  along  with  the 

contaminant  material  balance,  the  mass  balance  becomes 
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which  leads  to 
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Equations  (B-81)  and  (B-83)  are  combined  to  give 
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Furthermore,  Equation  (B-84)  Is  assumed  to  apply  when  (p  -  p^ 
constant . 


When  heat  transfer  from  the  surface  Is  present,  vertical  mixing 
will  be  enhanced  by  the  convection  turbulence  due  to  heat  transfer.  Zeman 
and  Tennekes  (Reference  B-16)  model  the  resulting  vertical  turbulent 
velocity  as 


w 


(B-85) 


where  w^  Is  the  convective  scale  velocity  described  as 
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If  u  Is  is  evaluated  at  H^ppt 
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where 
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and  is  the  temperature  obtained  from  the  energy  balance  of  Eqiiations 

(B-102)  and  (B-103).  Equation  (B-84)  is  modified  to  account  for  this 
enhanced  mixing  by 


[  ^l'^eff‘*eff  ] 
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where  Ri'  -  Ri.  — 

*  *  w 


Although  derived  for  two-dimensional  dispersion,  this  is  extended 
for  application  to  a  denser-than-air  gas  pltime  which  spreads  laterally  as  a 
density  intrusion: 
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where  the  plume  effective  half  width  is  defined  by 

B___  -  b  +  S 
EFF  2  y 

and  determined  using  the  gravity  intrusion  relation 


(B-91) 
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The  average  transport  velocity  in  the  plume  Is  defined  by 


and  the  lateral  spread  of  the  cloud  is  modeled  by 


(B-93) 


(B-94) 


2.  Horizontal  Dispersion 

The  crosswind  similarity  parameter  Sy(x)  is  also  determined  by 
requiring  that  it  satisfy  the  diffusion  equation 
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with  the  horizontal  turbulent  diffusivity  given  by 
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When  b  -  0,  S  •  Jl  a  ,  where  a  is  the  similarity  parameter  correlated 

y  y  y  a 

by  Pasquill  (Reference  B-16)  in  the  form  ••  6x  .  Furthermore,  Equations 
(B-95)  and  (B-96)  require  that 
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where  Equation  (B-98)  Is  also  assumed  applicable  for  determining  when  b 
Is  not  zero. 

At  the  downwind  distance  where  b  -  0,  the  crosswind  concentra¬ 
tion  profile  Is  asstimed  Gaussian  with  given  by 

Sy  -  yi  5(x  +  (B-99) 

where  x  Is  a  virtual  source  distance  determined  as 

V 

Sy(x^)  -  72  6<x^  +  x^)^  (B-lOO) 

The  gravity  spreading  calculation  Is  terminated  for  x  >  x^. 

For  a  steady  plume,  the  centerline  concentration  c  Is  determined 

c 

from  the  material  balance 


cu^dydz  -  2c^ 
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where  E  Is  the  plume  source  strength. 
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3 .  Energy  Balance 

For  some  simulations  of  cryogenic  gas  releases,  heat  transfer  to 
the  plume  In  the  downwind  dispersion  calculation  may  be  important,  particu¬ 
larly  in  low  wind  conditions.  The  source  calculation  determines  a  gas/alr 
mixture  initial  condition  for  the  downwind  dispersion  problem.  Air 
entrained  into  the  plume  is  assumed  to  mix  adlabatlcally .  Heat  transfer  to 
the  pl\une  downwind  of  the  source  adds  additional  heat.  This  added  heat  per 
unit  mass  Is  determined  by  an  energy  balance  on  a  uniform  cross  section 
as 


dx 
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where  q  is 
calculating 

dx 


determined  by  Equation  (B-37)  along  with  the  desired  method  of 
h^.  Equation  (B-102)  is  applied  when  b  -  0  and  is  extended  to 
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when  b  >  0.  Since  the  average  density  of  the  layer  p-^  cannot  be  determined 
until  the  temperature  (i.e.  Djj)  is  known,  a  trial  and  error  procedure  is 
required. 

Equations  (B-77),  (B-78),  (B-79).  (B-87) - (B-91) ,  (B-94),  and 
(B-98) - (B-103)  are  combined  with  an  equation  of  state  relating  cloud  den¬ 
sity  to  gas  concentration  and  temperature  and  are  solved  simultaneously  to 
predict  S^,  Sy,  c^,  and  b  as  functions  of  downwind  distance  beginning  at 
the  downwind  edge  of  the  gas  source. 


C.  CORRECTION  FOR  ALONG-WIND  DISPERSION 

Following  Colenbrander  (Reference  B-2),  an  adjustment  to  c^,  is  applied 
to  account  for  dispersion  parallel  to  the  wind  direction.  The  calculated 
centerline  concentration  Cj,(x)  is  considered  to  have  resulted  from  the 
release  of  successive  planar  puffs  of  gas  (Cj.(x)Ax)  without  any  dispersion 
in  the  x-dlrection.  If  it  is  assumed  that  each  puff  diffuses  in  the 
x-dlrectlon  as  the  puff  moves  downwind  Independently  of  any  other  puff  and 
that  the  dispersion  is  one -dimensional  and  Gaussian,  the  x- direction  con¬ 
centration  dependence  is  given  by 
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where  x^  denotes  the  position  of  the  puff  center  due  to  observer  i. 

After  Beals  (Reference  B-18),  the  x-direction  dispersion  coefficient 
is  assumed  to  be  a  function  of  distance  from  the  downwind  edge  of  the  gas 
source  (X  -  x  -  xg)  and  atmospheric  stability  given  by 

a  (X)  -  0.02  X^'^^  unstable,  x  >  130  m 

*  1.14 

-  0.04  X  neutral,  x  >  100  m  (B-105) 

“  0.17  X^’^^  stable,  x  >  50  m 

where  (X  -  x  -  x-)  and  a  are  in  meters.  The  concentration  at  x  is  then 

0  X 

determined  by  superposition,  i.e.,  the  contribution  to  at  a  given  x 

from  neighboring  puffs  is  added  to  give  an  x-direction  corrected  value  of 

c'.  For  N  observers, 
c 
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and  for  large  N, 
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The  corrected  centerline  concentration  c^  is  used  in  the  assumed 
profiles  in  place  of  c^.,  along  with  the  distribution  parameters  Sy,  S^, 
and  b. 


D.  DEGADIS  MODEL  INPUTS  AND  OUTPUTS 

As  implemented  under  VAX/VMS*,  DEGADIS  requires  three  areas  of  input 
information: 

•  simulation  definition 

•  numerical  parameters 

•  VAX/VMS  command  procedure  for  execution 

DEGADISIN  is  the  interactive  input  module  which  generates  the  simulation 
definition  from  user  responses.  An  example  input  session  is  included  in 
Section  D.3.  The  numerical  parameters  (convergence  criteria,  initial 
increments,  etc.)  are  supplied  to  DEGADIS  through  a  series  of  input  files. 
Although  these  numerical  parameters  are  easily  changed,  the  user  should 
need  to  change  these  only  rarely  with  the  exception  of  the  time  sort 
parameters .  The  VAX/VMS  command  procedure  used  to  execute  DEGADIS  is 
generated  in  DEGADISIN  by  default.  Additional  information  can  be  found  in 
Havens  and  Spicer  (Reference  B-1) . 

1 .  VAX/VMS  Command  Procedure 

The  VAX/VMS  command  procedure  generated  by  DEGADISIN  controls  the 
execution  of  images  for  the  simulation.  Image  execution  follows  one  of  two 
paths,  either  for  a  transient  release  or  for  a  steady-state  release. 


*VAX  and  VMS  are  registered  trademarks  of  Digital  Equipment  Corporation. 
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DEGADISIN  will  automatically  generate  the  appropriate  command  procedure. 
DEGADISIN  requires  a  simulation  name  be  specified.  The  simulation  name 
must  be  a  valid  VAX/VMS  file  name  without  a  file  eiiten^lon  and  is  desig¬ 
nated  herein  as  RUNNAME.  DEGADIS  will  use  this  file  name  with  standard 
extensions  for  input.  Interprocess  communication,  and  output. 

2.  Simulation  Definition 

DEGADISIN  is  an  Interactive  method  of  simulation  definition  where 
the  user  specifies  information  about  the  ambient  wind  field,  the  properties 
of  the  released  gas,  and  some  details  of  the  release.  A  summary  of 
required  input  information  is  included  in  Figure  B-5. 

The  ambient  wind  field  is  characterized  by  a  known  velocity  uq  at  a 
given  height  zq,  a  surface  roughness  and  the  Pasquill  stability  class. 
The  Pasquill  stability  class  along  with  the  desired  averaging  time  are  used 
to  estimate  values  of  the  lateral  similarity  parameter  coefficients  S  and 
(References  B-17  and  B-19  and  Appendix  D),  values  of  the  along-wlnd  simi¬ 
larity  coefficients  (Reference  B-18),  and  the  Monin-Obukhov  length  A  used 
by  Businger  et  al.  (Reference  B-20)  in  their  logarithmic  velocity  profile 
function.  The  Monin-Obukhov  length  is  then  used  to  calculate  the  friction 
velocity  u*.  Once  these  parameters  have  been  estimated  using  the  Pasquill 
stability  class,  the  user  has  the  option  of  interactively  changing  any  of 
these  to  better  describe  the  simulation.  In  addition  to  these  specifica¬ 
tions,  the  ambient  temperature,  pressure,  and  humidity  must  be  specified. 

The  properties  of  air  and  the  released  gas  are  used  to  evaliiate  the 
mixture  density  as  a  function  of  temperature  and  composition.  The  desired 
released  gas  properties  include  the  molecular  weight  MW^,  the  storage  tem¬ 
perature  (normal  boiling  point  for  cryogenic  gases)  Tq,  the  vapor  phase 
density  at  the  storage  temperature  and  ambient  pressure  pq,  and  two  con¬ 
stants  qj^  and  pj^  which  describe  the  heat  capacity  according  to  the  equation 


■ 

■ 

-1 

-  ,,.4 

Pi  Pi  ' 
T  -  Tq" 

C  (T)  -  (MW  ) 

Pc 

3.33  X  10  +  qj^ 

T  -  Tq 

where  C  (T)  is  the  mean  heat  capacity  (J/kg  K)  at  temperature  T.  Note 
Pc 

that  a  constant  heat  capacity  with  respect  to  temperature  can  be 
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Figure  B-5.  Summary  of  Simulation  Definition  Input  Information  for  DEGADIS. 


Including  Heat  Transfer 


Figure  B-5.  Summary  of  Simulation  Definition  Input  Information  for  DEGADIS 
(Continued) . 


Water  transfer  coefficient 


Summary  of  Simulation  Definition  Input  Information  for  DEGADIS 
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Figure  B-5.  Summary  of  Simulation  Definition  Input  Information  for  DEGADIS 
(Concluded) . 


obtained  by  setting  -  I.O  and  choosing  the  appropriate  value  for 
Representative  gas  properties  for  liquefied  natural  gas  (LNG)  as  methane, 
liquefied  petroleum  gas  (LFG)  as  propane,  pure  tinreacted  N02>  pure 
unreacted  N20^,  and  ammonia  (NH3)  are  Included  In  DEGADISIN.  Also  included 
are  the  lover  and  upper  flammability  limits  (LFL  and  UFL,  respectively), 
for  LNG  and  LPG. 

The  user  may  also  choose  to  calculate  the  mixture  density  as  a 
fxinctlon  of  composition  using  some  other  method.  This  mixture  density  is 
entered  In  the  program  as  If  the  release  were  Isothermal;  for  each  com¬ 
position,  the  program  requests  the  contaminant  mole  fraction,  the  contami¬ 
nant  concentration,  and  the  mixture  density.  For  ease  of  Input,  these 
values  may  be  entered  from  a  file  made  available  to  DEGACISIM. 

In  specifying  the  details  of  the  release,  the  user  must  choose  to 
simulate  the  release  as  transient  or  steady-state.  For  both  release  types, 
the  area  source  Is  assumed  circular.  The  source  radius  and  emission  rate 
must  be  specified  for  a  steady-state  release  only  once,  while  these  must  be 
specified  as  a  function  of  time  for  transient  releases  (either  Inter- 
I  actively  or  by  file).  For  transient  releases,  the  user  must  specify  the 


Initial  amount  of  gas  present  over  the  source  (in  order  to  simulate,  for 
example.  Instantaneous  releases  such  as  the  Thomey  Island  Trials). 

Figure  B-5  svimmarlzes  the  simulation  Information  gathered  by 
DEGADISIN  contained  In  the  RUNNAME  file  with  extension  INF.  The  structure 
of  RUNNAME. INF  is  Illustrated  In  Figure  B-6.  At  this  point,  RUNNAME. INF 
may  be  edited  to  correct  any  misinformation  entered  during  the  input 
session.  Note  that  care  must  be  exercised  when  editing  RUNNAME. INF  due  to 
the  fact  that  Information  contained  In  the  file  can  be  different  depending 
on  the  answered  questions  (e.g.  steady-state  versus  transient  simulation). 


Desert  Tortoise 

4.  Pressurized 

ammonia  release  by  LLNL. 

Released  mass  rate  of  100  kg/s; 

source  radius 

from  assumed  .H. 
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Figure  B-6.  DT4.1NP  Listing. 
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(The  reverse  of  this  page  is  blank.) 


3.  Example  Input  Sessions 

The  input  procedures  for  simulation  of  the  transient  release  and 
the  steady-state  release  are  very  similar.  Therefore,  only  the 
specification  of  a  steady- state  release  (DT4)  has  been  included.  In  the 
point  by-point  discussion  of  the  input  procedure,  note  the  following: 

(*)  The  file  name  specification  RUNNAME  must  satisfy  system 
restrictions . 

(*)  A  line  terminator  (normally  a  carriage  return)  must  end  every 
line  entered  by  the  user. 

(*)  When  DEGADISIN  requests  the  user  to  choose  an  option,  all 

acceptable  responses  are  a  single  character  (capital  or  lower 
case) .  The  default  responses  are  denoted  by  a  capital  letter 
inside  angle  brackets  (e.g.  <N>) .  When  applicable,  a  menu  of 
acceptable  responses  is  included  inside  parentheses. 

(*)  For  nvimerical  responses,  a  comma,  space,  tab,  or  line 
terminator  (carriage  return)  may  separate  the  numbers. 

(*)  When  a  file  is  used  as  input  (i.e.,  for  the  density  or 

transient  source  input),  DEGADISIN  reads  the  same  information 
from  the  file  which  would  be  entered  at  the  temninal  in  the 
same  order  and  in  the  same  format. 
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The  current  settings  of  pertinent  atmospheric  parameters  are 
displayed  in  this  list.  If  any  of  these  are  to  be  changed,  the  first 
letter  of  the  parameter  to  be  changed  is  entered.  The  value  of  the 
Monln-Obukhov  length  is  to  be  changed.  Note  that  the  default- - 
indicated  by  <N>--ls  No  for  no  changes. 
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(  No  ,  liel  t  a  ,  Beta  ,  Length , 


The  value  of  the  Monln-Obukhov  length  Is  indicated  to  be  infinite 
(which  is  the  default  value  for  a  Pasquill  D  stability) .  Any  value 
of  the  Monln-Obukhov  length  can  be  specified,  but  an  infinite  value 
is  a  special  case  which  can  be  obtained  for  any  Pasquill  stability 
class  by  setting  the  Monin-Obukhov  length  to  0.  The  desired  Monln- 
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Note;  For  infinity,  ML  ^  0.0 

Enter  the  desired  Monin-Obukhov  length;  (m)  45. 
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The  generated  INP  file  for  DT4  is  shown  in  Figure  B-6.  If  necessary, 
the  user  may  edit  the  INP  file  before  beginning  the  simulation. 

4.  Example  Simulation  Output 

After  proper  completion  of  the  model,  0T4.LIS  contains  the  output 
listing  for  the  steady-state  release.  A  point -by -point  discussion  of  the 
output  follows. 
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5.  Model  Limitations  and  Cautions 

DE6ADIS  model  application  should  be  limited  to  the  description  of 
atmospheric  dispersion  of  denser -than -air  gas  releases  at  ground  level  onto 
flat,  unobstructed  terrain  or  water.  Application  to  releases  from  sources 
above  ground  level  (e.g.  overflow  from  dikes)  would  be  expected  to  give 
conservative  predictions  of  the  downwind  hazard  zones ,  but  this  has  not 
been  verified. 

The  dispersion  of  a  denser -than -air  gas  by  the  action  of  the  wind 
assumes  the  maintenance  of  a  wind  velocity  profile  in  the  gas  cloud  or 
plume  whose  characteristics  are  determined  by  the  approach  wind  flow 
(upwind  of  the  release) .  The  treatment  of  vertical  momentum  transfer 
invokes  the  assumption  of  of  a  logarithmic  vertical  velocity  profile,  which 
is  in  turn  curve-fitted  to  a  power  law  vertical  velocity  profile.  DEGADIS 
also  uses  similarity  forms  for  the  vertical  profile  of  gas  concentration  in 
the  cloud,  and  the  vertical  profile  is  dependent  on  the  power  law  exponent 
a  used  in  the  representation  of  the  velocity  profile.  The  vertical 
velocity  profile,  which  is  directly  related  to  the  air  entrainment  velocity 
into  the  cloud,  is  dependent  on  the  factors  which  determine  the  structure 
of  the  atmospheric  boundary  surface  layer,  wind  speed,  surface  roughness, 
and  atmospheric  stability.  Consequently,  the  representations  of  the 
vertical  velocity  and  concentration  profiles  in  DEGADIS  are  subject  to 
similar  limitations  as  in  other  descriptions  of  the  surface  layer.  Table 
B-2  indicates  typical  recommended  surface  roughness  values.  Table  B-1 
indicates  logarithmic  wind  velocity  profile  corrections  for  different 
atmospheric  stabilities,  along  with  typical  values  of  the  wind  profile 
power  law  exponent  a  determined  in  DEGADIS. 

Demonstration  of  the  model  has  been  primarily  directed  to  the 
prediction  of  hazard  extent  defined  by  gas  concentrations  in  the 
hydrocarbon  flammable  limit  range  (~1  to  20  percent) .  Based  on  the 
simulations  of  field  experiments  presented  in  Havens  and  Spicer  (Reference 
B-1),  the  ratio  of  observed  distance  to  calculated  distance  for  a  given 
time-averaged  concentration  level  (OBS/PRE)  ranged  from  0.73  to  0.96  for 
the  5  percent  level  9  out  of  10  times  (i.e.,  90  percent  confidence 
interval).  For  the  2.5  percent  level,  (OBS/PRE)  ranged  from  0.82  to  1.03 
for  a  90  percent  confidence  interval.  For  the  1  percent  level,  (OBS/PRE) 
ranged  from  0.95  to  1.24  for  a  90  percent  confidence  Interval.  If  for  a 


TABLE  B-2.  REPRESENTATIVE  VALUES  OF  SURFACE  ROUGHNESS  FOR  A  UNIFORM 

DISTRIBUTION  OF  THESE  TYPES  OF  GROUND  COVER  (REFERENCE  B-21) 


Surface  Roughness 
(m) 

Height  of 

Gro\ind  Cover  (ra) 

Ice 

0.00001 

Smooth  mud  flats 

0.00001 

Smooth  snow  on  short  grass 

0.00005 

Snow 

0.00005  to  0.0001 

Sand 

0.0003 

Smooch  desert 

0 . 0003 

Snow  surface,  natural 

prairie 

0.001 

Soils 

0.001  to  0.01 

Short  grass 

0.003  to  0.01 

0.02  to  0.1 

Mown  grass 

0.002 

0.015 

0.007 

0.03 

Long  grass 

0.04  to  0.10 

0.25  to  1. 

Agricultural  crops 

0.04  to  0.20 

-0.40  to  2 . 

Orchards 

0.50  to  1. 

—5 .  to  10 . 

Deciduous  forests 

1 .  to  6 . 

—10 .  to  60 . 

Coniferous  forests 

1 .  to  6 . 

-10.  to  60. 

given  release  scenario  Che  calculated  distance  to  the  2.5  percent  average 
concentration  level  was  120  meters,  the  distance  to  the  2.5  percent  average 
concentration  for  9  out  of  10  realizations  of  the  same  release  would  be 
expected  to  range  between  98  meters  and  124  meters. 
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APPENDIX  C 

DISCUSSION  OF  THE  SENSITIVITY  OP  DEGADIS  TO  UNCERTAINTY  IN 
INPUT  PARAMETERS  NOT  ROUTINELY  AVAILABLE  AT 
OPERATIONAL  USAF  SITES 


The  purpose  of  this  appendix  is  to  point  out  the  effect  of  uncertain¬ 
ties  in  DEGADIS  input  parameters  on  the  output  of  DEGADIS.  The  scope  of 
this  investigation  is  aimed  at  input  parameters  which  are  not  normally 
available  at  operational  USAF  sites.  For  this  analysis,  the  following  data 
are  assumed  to  be  available: 

•  windspeed  at  a  given  elevation  (assumed  to  be  near  10  meters) 

•  ambient  temperature,  pressure,  and  humidity 

(The  ambient  humidity  is  assumed  to  be  available  from  the  dewpoint 
temperature.)  In  addition  to  the  above  information,  DEGADIS  requires 
further  information  to  define  a  simulation  as  follows; 

•  Pasquill  stability  category 

•  surface  roughness 

•  initial  contaminant  density 

•  contaminant  release  rate 

The  release  conditions  of  Desert  Tortoise  4  (DT4)  will  be  used  to 
illustrate  the  effect  of  tmcertainties  in  each  of  these.  For  comparison 
purposes,  the  DT4  conditions  are  the  same  as  those  In  the  body  of  this 
report  with  two  exceptions.  The  first  is  the  Monln-Obukhov  length. 

Because  LLNL  measured  the  on-site  velocity  and  temperature  profiles,  an 
estimate  of  the  Monin-Obukhov  length  was  available.  This  estimate  was  used 
in  the  simulations  for  comparison  with  the  available  data.  Since  the 
Monin-Obukhov  length  is  not  a  parameter  which  is  normally  available  and 
since  DEGADIS  can  estimate  the  Monin-Obukhov  length  from  the  Pasquill 
stability  category  and  the  surface  roughness,  the  Monin-Obukhov  length  will 
be  estimated  by  DEGADIS  (as  outlined  in  Table  B-1).  The  second  is  the 
mixture  density  specification  as  a  function  of  contaminant  concentration. 
For  releases  where  heat  transfer  is  not  important,  assuming  a  linear 
relationship  between  contaminant  concentration  (in  kg/m-’)  and  mixture 
density  may  be  adequate  for  hazard  assessment  as  shown  by  comparison  of 
DEGADIS  predictions  using  the  conditions  of  DT4  (Section  III) .  (If  heat 
transfer  is  important,  the  effect  of  heat  transfer  would  be  expected  to 
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enhance  dispersion  so  that  a  simulation  vhlch  Ignores  heat  transfer  would 
be  expected  to  be  conservative.)  For  all  comparisons,  a  linear 
relationship  Is  assumed  between  contaminant  concentration  and  mixture 
density. 

The  Pasqulll  stability  category  and  site  surface  roughness  are  needed 
in  addition  to  the  normally  available  Information  outlined  above  to 
complete  the  description  of  the  ambient  meteorological  conditions  of  the 
release.  The  Pasqulll  stability  and  surface  roughness  are  used  to  estimate 
the  power  law  wind  profile  parameters  (a),  the  Monln-Obukhov  length  (A), 
and  the  coefficients  for  the  lateral  and  along-wlnd  dispersion  parameters 
(ay  and  .  The  parameter  \  determines  the  shape  of  the  ambient  velocity 
profile  (Equation  B-73)  and  Is  used  to  estimate  and  a.  The  parameter  a 
Is  used  In  the  vertical  concentration  profile  and  the  vertical  velocity 
profile.  The  Oy  coefficients  are  used  In  the  determination  of  Syi  the 
coefficients  are  only  used  In  transient  simulations  for  the  x-dlrectlon 
dispersion  correction. 

There  are  several  methods  vhlch  can  be  used  to  estimate  the  Pasqulll 
stability  category.  An  estimate  of  the  Pasqulll  stability  category  from 
the  cloud  cover  and  wind  speed  can  be  obtained  as  shown  In  Table  C-1 
(Reference  C-1  or  C-2) .  If  Information  on  the  horizontal  wind  direction  Is 
available,  a^  can  be  estimated  from  the  maximum  deviation  In  the  horizontal 
wind  direction  (3a^  maximum  deviation);  the  corresponding  stability 
categories  are  shown  In  Table  C-2  (Reference  C-1  or  C-2). 

A  comparison  of  DT4  simulations  using  D,  E,  and  F  stabilities  is  shown 
in  Figure  C-1.  A  shown,  the  DEGADIS -predicted  concentrations  at  a  given 
distance  for  each  stability  class  are  within  a  factor  of  two  for  concen¬ 
trations  greater  than  about  1000  ppm. 

The  site  surface  roughness  Is  normally  determined  by  analysis  of  the 
velocity  profile  measured  under  a  given  set  of  circumstances  for  a  specific 
site.  Estimates  of  the  surface  roughness  for  a  site  are  normally  con¬ 
sidered  to  be  only  a  function  of  the  ground  cover;  representative  values 
for  homogeneous  surfaces  are  shown  in  Table  C-3.  As  a  first  approximation, 
the  surface  roughness  can  be  estimated  as  one-fifth  to  one-tenth  of  the 
terrain  height  for  homogeneous  terrain.  For  urban  and  suburban  areas, 
Lettau  (Reference  C-3)  suggests  Zj^  >  hA  /2A  for  a  uniform  distribution  of 
buildings  where  h  Is  the  building  height,  A  Is  the  area  normal  to  the  mean 


TABLE  C-1.  PASQUILL  STABILITY  CATEGORY  AS  A  FUNCTION  OF  tfINDSPEED  AND 
CLOUD  COVER  (REFERENCE  C-1  AND  REFERENCE  C-2) 


Surface 

Vindspeed 

(m/s) 

Day t ime  Insolation _ 

Nittht 

Strong* 

Moderate^ 

Slight‘d 

Thin  Overcast 

or 

>4/8  Low  Cloud 

<3/8 

Cloud 

<2 

A 

A-B 

B 

t 

CM 

A-B 

B 

C 

E 

F 

3-5 

B 

B-C 

C 

D 

E 

5-6 

C 

C-D 

D 

D 

D 

>6 

C 

D 

D 

D 

D 

S-BOO  V/m^ 

^600  to  300  W/m2 
*^<300  W/m^ 


TABLE  C-2.  PASQUILL  STABILITY  CATEGORY  AS 
A  FUNCTION  OF  HORIZONTAL  WIND 
DIRECTION  FLUCTUATION  (FROM 
REFERENCE  C-1  AND  REFERENCE  C-2) 


Pasquill  Stability  Category 

og  (at  10  m) 

A 

25 

B 

20 

C 

15 

D 

10 

E 

5 

F 

2.5 

10^' 


TABLE  C-3.  REPRESENTATIVE  VALUES  OF  SURFACE  ROUGHNESS  FOR  A  UNIFORM 

DISTRIBUTION  OF  THESE  TYPES  OF  GROUND  COVER  (REFERENCE  C-4) 


Surface  Roughness 
(m) 

Height  of 

Ground  Cover  (m) 

Ice 

0.00001 

Smooth  mud  flats 

0.00001 

Smooth  snow  on  short  grass 

0.00005 

Snow 

0.00005  to  0.0001 

Sand 

0.0003 

Smooth  desert 

0.0003 

Snow  surface ,  natural 

prairie 

O.OOl 

Soils 

0.001  to  0.01 

Short  grass 

0.003  to  0,01 

0.02  to  0.1 

Mown  grass 

0.002 

0.015 

0.007 

0.03 

Long  grass 

0.04  to  0.10 

0.25  to  1. 

Agricultural  crops 

0,04  to  0.20 

-0.40  to  2 . 

Orchards 

0.50  to  1. 

-5.  to  10. 

Deciduous  forests 

1 .  to  6 . 

-10.  to  60. 

Coniferous  forests 

1 .  to  6 . 

-10.  to  60. 

wind  direction,  and  a'  Is  the  building  lot  area.  Of  course,  the  area 
downwind  of  the  source  may  be  characterized  by  distinct  terrain  types 
having  different  surface  roughnesses;  for  such  conditions,  separate  calcu¬ 
lations  using  the  minimum  and  maximum  roughnesses  will  act  as  upper  and 
lower  bounds  for  the  simulation.  It  should  be  noted  that  DEGADIS  Is  based 
on  the  asstunptlon  of  dispersion  on  unobstructed,  flat  terrain.  When  the 
surface  roughness  used  In  DEGADIS  becomes  a  significant  fraction  of  the 
depth  of  the  dispersing  layer,  the  underlying  assumptions  In  DEGADIS  may  no 
longer  be  satisfied. 
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A  comparison  of  DT4  simulations  using  surface  roughnesses  of  0.3,  0.03, 
and  0.003  meters  is  shown  in  Figure  C-2.  As  shorn  the  DEGADIS -predicted 
concentrations  at  a  given  distance  are  reduced  by  a  factor  of  about  two  for 
each  order  of  magnitude  increase  in  surface  roughness. 

For  releases  of  materials  which  are  not  significantly  colder  than  the 
ambient  temperature  (say,  greater  than  about  200  K) ,  the  most  important 
property  of  the  released  gas  is  the  initial  gas  density  and  the  density  of 
the  contaminant/air  mixture  as  a  function  of  contaminant  concentration. 

For  releases  where  heat  transfer  is  not  important,  assuming  a  linear 
relationship  between  contaminant  concentration  (in  kg/m^)  and  mixture 
density  may  be  adequate  for  hazard  assessment.  (If  heat  transfer  is 
important,  the  effect  of  heat  transfer  would  be  expected  to  enhance 
dispersion  so  that  a  simulation  which  ignores  heat  transfer  would  be 
expected  to  be  conservative.) 

A  comparison  of  DT4  simulations  using  different  initial  contaminant 
densities  ranging  from  20  percent  below  (2.98  kg/m^)  to  20  percent  above 
(4.47  kg/m^)  the  original  initial  contaminant  density  (3.72  kg/m^)  showed 
that  the  DEGADIS-predicted  concentration  as  a  function  of  distance  was 
essentially  unchanged.  A  0T4  simulation  using  an  initial  contaminant 
density  of  1.86  kg/m^  showed  a  variation  of  about  20  percent  in  the 
DEGADIS-predicted  maximum  concentration  at  a  given  distance.  Note  that 
this  density  is  still  much  greater  than  the  ambient  air  density  (1.02 
kg/m^) . 

The  last  category  to  be  examined  is  the  release  method  of  the  material. 
For  ground-level  releases,  two  considerations  are  important  including  the 
rate  of  release  and  the  diameter  of  release.  As  pointed  out  herein,  the  DT 
simulations  were  found  to  be  insensitive  to  the  choice  of  the  source 
diameter.  Furthermore,  the  effect  of  the  source  rate  can  be  seen  by 
comparing  the  simulations  for  each  test  for  the  two  mass  rates  simulated 
(the  released  mass  rate  and  the  mass  rate  passing  the  800-meter  sensor 
array).  As  shown  in  these  steady-state  simulations,  the  concentration  at  a 
given  distance  for  a  given  set  of  release  conditions  is  proportional  to  the 
contaminant  evolution  rate.  This  is  the  same  behavior  shown  by  the  steady- 
state  Gaussian  plume  model. 


10^ 


Maximum  Concentration  (%) 


0.003  m 
0.03  m 


0.3  m 


Figure  C-2.  Comparison  of  DEGADIS-Predicted  Maximum 
Concentration  as  a  Function  of  Distance 
for  Surface  Roughnesses  of  0.3,  0.03, 
and  0.003  m  for  the  Test  Case. 
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From  the  above  analysis,  the  Inputs  which  have  the  largest  effect  on 
the  DEGAOIS -predicted  distance  to  a  given  concentration  level  are  the 
source  evolution  rate,  sturface  roughness,  and  the  Pasqulll  stability  class 
for  steady-state  simulations.  In  addition,  the  Initial  contaminant  density 
appears  to  be  Insignificant  as  long  as  It  Is  above  some  threshold  value. 
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APPENDIX  D 

USING  DIFFERENT  TIME-AVERAGING  PERIODS  IN  DEGADIS 


Because  of  the  nature  of  the  exposure  of  Individuals  to  airborne  toxic 
materials,  It  is  necessary  to  be  able  to  predict  the  exposure  to  a  given 
compound  over  a  specified  time  period.  Estimation  of  the  maximum  downwind 
concentration  for  a  given  averaging  time  using  DEGADIS  has  been  addressed 
by  Spicer  (Reference  D-l);  the  following  Is  a  summary  of  that  work. 

Although  the  effect  of  averaging  time  on  the  maximum  downwind  concentration 
for  steady  releases  Is  still  open  to  some  question,  the  most  Important 
effect  will  be  assumed  to  be  the  result  of  plume  meander  for  modeling 
purposes.  For  transient  releases,  the  most  Important  effect  will  be 
assumed  to  be  the  duration  of  the  release  In  addition  to  the  mechanism  of 
plume  meander.  (If,  for  example,  the  duration  of  the  release  Is  much 
shorter  than  the  averaging  time,  then  the  time-averaged  concentration  will 
be  reduced.)  In  the  following,  the  Implications  of  the  averaging  time  on 
the  Gaussian  plume  model  are  discussed,  and  a  way  of  using  different 
averaging  times  In  DEGADIS  for  risk  assessment  purposes  Is  presented. 

The  effect  of  averaging  time  on  the  maximum  downwind  concentration  for 
steady  releases  Is  generally  recognized  as  being  a  result  of  plume  meander. 
The  plume  from  a  steady  release  of  a  passive  gas  would  be  expected  to  move 
downwind  and  meander  with  the  ambient  wind.  At  any  point  In  time,  the 
centerline  of  the  plume  would  not  necessarily  correspond  with  the  mean  wind 
direction.  After  averaging  the  plume  boundaries  over  some  time  period,  the 
centerline  of  the  plume  should  more  closely  correspond  with  the  mean  wind 
direction.  This  behavior  Is  reflected  In  the  values  of  the  dispersion 
parameters  used  In  the  Gaussian  plume  model;  Turner  (Reference  D-2)  and 
Beals  (Reference  D-3)  report  different  dispersion  parameters  for  10-mlnute 
average  plume  behavior  and  Instantaneous  or  puff  behavior.  (It  should  be 
noted  that  estimates  of  puff  coefficients  are  usually  derived  from  obser¬ 
vations  of  continuous  releases  which  are  analyzed  as  though  the  plume  Is  a 
continuous  train  of  superimposed  puffs  (Reference  D-4).)  Because  of  plume 
meander,  the  lateral  Gaussian  dispersion  parameter  Oy  depends  on  the 
averaging  time,  while  the  value  of  Is  essentially  unaffected  by  the 
averaging  time  (Reference  D-2) .  Several  Investigators  have  also  looked  at 
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the  effect  of  avereging  time  by  examining  the  ratio  of  maxlimim 
concentrations  for  different  averaging  times  for  continuous  releases;  the 
usual  form  of  the  function  specified  is 

where  c,^^y.^(x;t-| )  is  the  maximum  concentration  associated  with  averaging 
time  t]^,  Cn,«.g(a;t^)  is  the  maximum  concentration  associated  with  averaging 
time  t2.  and  p  is  some  power.  Turner  (Reference  D-2)  proposed  that  0.17  < 
p  ^  0.20  based  on  reports  by  Stewart.  Gale,  and  Crooks  (Reference  D-5)  and 
Cramer  (Reference  D-6)  (among  others)  for  averaging  times  from  about  3 
seconds  to  about  half  an  hour.  Herman  (Reference  D-7)  estimated  p  to  be 
0.5  +  0.2  for  averaging  times  from  1  hour  to  1  year.  These  findings  are  in 
essential  agreement  with  Hino  (Reference  D-8)  who  found  Chat  p  -  0.2  for 
averaging  times  less  than  10  minutes,  and  p  >  0.5  for  times  greater  than  10 
minutes  (with  some  dependence  on  the  atmospheric  stability).  Of  course,  if 
the  value  of  is  independent  of  averaging  time  as  proposed  by  Beals 
(Reference  D-3),  Equation  (D-1)  can  be  restated  as 

<^y<*:tj^)/<7y(x;t2)  - 

where  ory(x;tj^)  and  <jy(x;t2)  are  the  values  of  associated  with  the 
averaging  times  tj^  and  t2,  respectively.  Obviously,  Equation  (D-2)  will 
not  be  appropriate  if  the  averaging  time  is  taken  to  be  zero  (i.e.,  for  a 
puff).  However,  Eq\iation  (D-2)  can  be  solved  for  an  effective  averaging 
time  associated  with  a  puff  value  of  Wy.  Using  the  values  reported  by 
Turner  (Reference  D-2),  if  the  puff  and  10-minute  pliime  values  for  for  D 
stability  are  used  in  Equation  (D-2),  the  effective  averaging  time  for  the 
puff  coefficient  is  about  20  seconds.  This  result  indicates  that  the  width 
of  a  steady-state  plume  would  not  be  expected  to  vary  significantly  due  to 
meander  over  any  given  -20  second  period.  (Note  that  the  last  statement 
would  not  be  approporiate  if  made  about  the  maximum  concentration.) 

Puff  values  of  Oy  reported  by  Turner  (Reference  D-2)  and  Beals 
(Reference  D-3)  are  given  for  the  general  stability  categories  stable, 
neutral,  and  unstable.  If  all  of  the  puff  Oj  values  are  assigned  the  same 
minimum  averaging  time  of  20  seconds,  the  puff  ay  values  for  the  stable 
category  approximately  correspond  to  F  stability,  and  the  puff  ay  values 
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for  the  unsteble  category  approxiffletely  correspond  to  B  stability,  if  the 
value  of  ffy  la  parameterized  as 


where  <7y  and  z  are  in  meters,  then  the  value  of  can  be  approximated  as 
being  the  same  for  the  plume  and  puff  values.  (Seinfeld  (Reference  D-A) 
used  -  0.894  for  the  plume  Oy  values.)  Using  the  same  power  /9,  the 
parameterlzatlons  for  plume  and  puff  Oy  values  are  shown  In  Table  D-1. 


TABLE  D-1.  COEFFICIENT  S  IN  GAUSSIAN  DISPERSION  MODEL  FOR  USE  IN 
CTy  -  fix^  WITH  0  -  0.894  AND  <7y  AND  X  IN  METERS 


Stability  Class 

Averaging 

Time 

A 

B 

C 

0 

E 

F 

10  min 

0.443 

0.324 

0.216 

0.141 

0.105 

0.071 

20  s  or  less 

0.224 

0.164 

0.109 

0.071 

0.053 

0.036 

In  the  steady  downwind  dispersion  phase  of  DEGADIS,  Equation  (D-3)  is 
Implemented  to  detenalne  the  lateral  dispersion  parameter  Sy  as  a  function 
of  distance.  (When  the  central  horlzontally-homogeneous  section  disappears 
(b  -  0),  ‘^y  "  Sy/y2.)  Due  to  the  assumed  profile  of  the  area  release  In 
DEGADIS,  the  Initial  value  of  Sy  Is  zero.  Downwind  of  the  source,  the 
gtowth  of  Sy  satisfies  the  diffusion  eqtiatlon  and  a  lateral  dispersion 
coefficient  consistent  with  the  specification  of  Oy  given  by  Eqtiatlon 
(D-3).  The  rate  of  growth  of  the  effective  plume  half -width  Bg££  Is 
determined  by  lateral  gravity  spreading  when  b  Is  nonzero.  (The  value  of  b 
is  actually  determined  by  b  -  B^f^  -  Sy/2.)  At  some  distance  downwind 

of  the  source  (x,.) ,  b  goes  to  zero,  then,  B^ff  -  Tx  Sy/2;  the  effective 
plume  width  Is,  thereafter,  determined  using  Equation  (D-3),  the  value  of 


Sy(x^),  and  a  virtual  source  distance  (x^)  determined  from  Sy.(x^)  -  J2  S 
(x^  -  x^)^-  Since  the  lateral  gravity  spreading  Is  treated  as  being 
Independent  of  the  plume  meander,  the  only  consideration  of  averaging  time 
for  the  steady  DEGADIS  downwind  dispersion  phase  Is  the  appropriate  choice 
of  the  constants  to  be  used  In  Equation  (D-3).  In  DEGADIS,  the  default 
value  for  Is  0.894;  the  default  value  for  S  Is  chosen  based  on  an 
averaging  time  of  20  seconds  or  less. 

For  transient  releases,  DEGADIS  uses  an  observer  scheme  to  analyze  the 
release.  (A  release  Is  considered  transient  when  the  release  duration  Is 
less  than  the  time  required  for  the  gas  to  travel  to  the  position  of 
Interest.  When  considering  the  effect  of  averaging  time,  a  release  Is  also 
considered  transient  when  the  release  duration  Is  less  than  the  averaging 
time.)  A  series  of  observers  which  travel  with  the  (approximate)  mean 
advectlon  velocity  of  the  gas  are  released  upwind  of  the  source  and  travel 
downwind.  As  each  observer  travels  over  the  source,  time  averages  of  the 
pertinent  source  parameters  are  determined.  Vlth  this  averaged  source,  a 
steady  downwind  dispersion  phase  calculation  Is  performed  for  each 
observer.  Since  the  (approximate)  observer  velocity  Is  known,  the  observer 
position  as  a  function  of  time  can  then  be  determined,  and  then,  concentra¬ 
tion  as  a  function  of  distance  and  time  can  be  determined.  For  transient 
releases.  It  Is  necessary  to  average  the  predicted  concentration  over  time 
with  the  Integral 


c(x,y,z,t;t^^) 


1 

t 

av 


r' 

c(x,y,z,t) 


t-t 


av 


dt 


<D-4) 


where  t^.^  Is  the  averaging  time. 

The  following  methodology  Is  recommended  for  using  different  averaging 
times  In  DEGADIS: 

(1)  To  account  for  the  effect  of  plume  meander  on  time-averaged 

concentration,  the  appropriate  value  of  $  from  Table  D-1  should  be 
used  In  the  simulation  whether  the  release  can  be  considered  steady 
or  not. 
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(2)  If  the  release  cannot  be  considered  steady  due  to  either  considera¬ 
tion  outlined  above,  then  the  release  must  be  simulated  as  a 
transient  release  and  the  resulting  concentrations  averaged  as  in 
Equation  (D-4) . 


114 


REFERENCES  FOR  APPENDIX  D 


D-1.  Spicer,  T.O.,  Halng  Different  Time  Averaging  Periods  in  DEGADIS. 
Report  to  the  Exxon  Education  Foundation,  June  1987. 

D-2.  Turner,  D.B.,  Workbook  of  Atmospheric  Dispersion  Estimates.  USEPA 

999-AP-26,  D.S.  Environmental  Protection  Agency,  Washington  DC,  1970. 

D-3.  Beals,  G.A.,  A  Guide  to  Local  Dispersion  of  Air  Pollutants.  Air 
Weather  Service  Technical  Report  214,  April  1971. 

D-4.  Seinfeld,  J.H. ,  "Atmospheric  Diffusion  Theory,"  Advances  in  Chemical 
Engineering.  12,  1983,  pp.  209-299. 

D-5.  Stewart,  N.G. ,  Gale,  H.J.,  and  Crooks.  R.  N. ,  "The  Atmospheric 

Diffusion  of  Gases  Discharged  from  the  Chimney  of  the  Howell  Reactor 
BEPO,"  Int.  J.  Air  Poll..  1.  1958,  pp.  87-102. 

D-6.  Cramer,  H.E.,  "Engineering  Estimates  of  Atmospheric  Dispersal 
Capacity,"  Amer.  Ind.  Hvg.  Assoc.  J. .  20.  1959,  pp.  183-189. 

D-7.  Herman,  M.N..  "Estimating  Long-Term  Ground  Level  Concentrations  of 
SO2  from  Short-Term  Peak  Data,"  JAPCA .  30.  1980,  pp.  676-678. 

D-8.  Hino,  Mikio,  "Maximum  Gro\ind-Level  Concentration  and  Sampling  Time," 
Atmospheric  Environment.  2,  1968,  pp.  149-165. 


115 

(The  reverse  of  this  page  is  blank.) 


